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ABSTRACT 


This  report  is  concerned  with  the  application  of  frac¬ 
ture  mechanics  theory  to  the  problem  of  projectile  pene¬ 
tration  and  the  situation  where  an  excessive  amount  of 
energy  is  available  to  damage  the  target  material  in  the 
form  of  fracture  followed  by  crack  propagation.  The  con¬ 
ditions  under  which  this  phenomenom  takes  place  depend  on 
several  parameters  such  as  the  speed  and  shape  of  the  pro¬ 
jectile  and  the  geometry  and  material  of  the  target.  These 
conditions  have  been  simulated  in  the  laboratory.  The 
present  investigation  deals  only  with  the  development  of  an 
analytical  fracture  model  of  the  projectile  penetration 
problem  using  the  theory  of  continuum  mechanics. 


The  amount  of  energy  required  to  produce  a  crack  of 
certain  size  in  a  pre-stretched  plate  within  a  very  short 
period  of  time  is  estimated  by  solving  a  system  of  dual 
integral  equations.  Computer  programs  are  being  developed. 
In  addition,  the  energy  supplied  to  maintain  a  running 
crack  at  a  certain  speed  is  obtained  with  results  illus¬ 
trated  by  graphs.  The  effects  of  plasticity  and  rate 
sensitivity  of  the  material  are  also  discussed  in  this 
report . 
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SECTION  I 


INTRODUCTION 

One  class  of  problems  that  has  received  little  atten¬ 
tion  relative  to  its  potential  practical  importance  is  that 
of  the  dynamic  fracture  of  structures  due  to  impact  loadings. 
This  is  mainly  because  of  the  extreme  complexities  of  the 
problems  which  involve,  in  addition  to  the  inertia  of  the 
structural  system,  such  effects  as  plastic  deformation,  rate 
dependent  material  property,  non-linear  behavior,  etc. 

Lacking  in  particular  is  an  understanding  of  the  failure  of 
structural  members  damaged  by  projectile  penetration. 

From  a  thermodynamic  viewpoint,  the  energy  input  to  a 
target  material  is  divided  into  elastic  strain  energy,  plas¬ 
tic  energy  dissipated  in  the  form  of  heat,  kinetic  energy, 
and  the  energy  carried  by  the  projectile  leaving  the  target 
material.  However,  if  the  input  work  exceeds  the  dissipation, 
then  the  unbalanced  portion  of  the  energy  will  be  absorbed 
in  creating  new  surfaces  starting  at  the  edge  of  a  stress 
raiser  of  the  damaged  target.  This  portion  of  the  total 
energy  is  referred  to  in  fracture  mechanics  as  the  "surface 
energy"  which  is  also  associated  with  quantities  commonly 
known  as  the  energy  release  rate  or  crack  driving  force  6 
and,  indirectly,  the  stress-intensity  factor  k  The 

relationship  between  G  and  k  has  been  recently  established 

(3) 

by  Sih  for  all  the  three  fundamental  modes  of  dynamic 


fracture.  Hence,  a  knowledge  of  G  or  k  will  provide  the 
basis  for  a  more  efficient  assessment  of  how  cracks,  once 
initiated,  might  be  stopped,  and  provide  the  first  step 
toward  the  possibility  of  defining  a  dynamic  fracture  tough¬ 
ness  value  of  a  material  which,  in  essence,  determi nes  the 
criterion  for  arresting  running  cracks. 

The  overall  problem  of  projectile  penetration  is  con¬ 
siderably  more  complicated  since  it  involves  not  only  an 
estimate  of  G  or  k  but  the  proportions  and  transfer  rates 
among  all  the  energies  described  earlier  throughout  the 
impact  process.  Aside  from  the  fracture  energy,  it  is  also 
important  to  know  the  portion  of  the  total  energy  that  is 
dissipated  through  transformation  of  the  target  material 
which  undergoes  a  transition  from  the  fluid  to  solid  state 
during  and  after  impact.  A  study  of  the  material  flow  pro¬ 
cess  will  involve  the  assumption  of  an  appropriate  consti¬ 
tutive  equation  and  a  numerical  scheme  which  necessitates 
large  scale  computations.  The  final  configuration  of  the 
damaged  target,  such  as  hole  diameter  a .. d  the  extent  or  speed 
of  petalling  cracks,  depends  on  the  distribution  of  the  re¬ 
maining  input  energy. 

The  material  flow  process  concerned  with  the  formation 
of  holes  in  plates  or  the  perforation  of  impacted  plates 
has  been  studied  by  Taylor  Chen  and  others 

This  report  deals  mainly  with  crack  formation  under  impact 
t  and  the  conditions  under  which  the  cracks  or  flaws  might 
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trigger  catastrophic  failure  of  the  target  material.  To 
this  end,  the  first  year  effort  has  been  devoted  to  assess 
the  amount  of  energy  dissipated  by  the  sudden  appearance  of 
a  crack  in  a  prestretched  target  material,  a  crack  running 
at  a  constant  speed,  and  the  effect  of  plastic  deformation 
ahead  of  a  moving  crack.  Other  considerations,  such  as  crack 
acceleration  and  rate  sensivity  of  the  materi al,  wi 1 1  be 
treated  in  the  future. 

A  general  expression  for  the  dynamic  energy  release 
rate  of  a  crack  undergoing  arbitrary  motion  is  developed  in 
Section  II.  This  is  essentially  the  work  done  for  opening 
a  line  segment  ahead  of  the  crack.  For  constant  crack  ve¬ 
locity,  the  expression  reduces  to  that  of  the  path  indepen- 

(3) 

dent  integral  derived  in  . 

The  problem  of  transient  loading  on  a  crack  is  considered 
in  Section  III.  Use  is  made  of  integral  transforms  coupled 
with  the  technique  of  Cagniard  v  1  used  for  solving  geo¬ 
physics  problems.  There  are  several  ways  that  this 
model  may  be  of  use.  For  example,  perforation  of  a  stressed 
plate  (or  pressure  vessel)  by  spall  debris  could  result  in 
what  may  be  approximated  as  the  sudden  appearance  of  a  crac'- 
and  subsequently  lead  to  catastrophic  failure.  In  an  un¬ 
stressed  plate  the  projectile  itself  could  create  and,  over 

a  short  period  of  time,  load  a  crack.  Such  a  crack  can  be 

( g ) 

seen  in  Figure  78  of  Reference  v  .  The  results  in  terms 
of  dynamic  stress-intensity  factors  are  plotted  as  a  function 


of  time. 
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In  Section  IV,  several  models  of  propagating  cracks 
are  considered.  Through  use  of  a  simple  model,  in  which  the 
plastic  yielding  is  assumed  to  occur  in  a  narrow  region 
ahead  of  the  running  crack,  the  amount  of  energy  dissipated 
in  the  various  running  crack  models  is  calculated.  This 
information  is  essential  in  estimating  the  fracture  energy 
of  the  projectile  penetration  process. 

Thus,  the  overall  objective  of  understanding  the  pro¬ 
jectile  penetration  problem  will  be  accomplished  by  first 
solving  a  number  of  idealized  problems  where  the  various 
effects  can  be  treated  separately.  Evaluation  of  the  pro¬ 
blem  as  a  whole  requires  a  coordinated  effort  of  research 
in  several  areas  which  are  still  being  investigated. 
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SECTION  II 

DYNAMIC  CRACK  DRIVING  FORCE 

An  expression  is  derived  here  for  the  energy  release 
rate  associated  with  a  moving  crack  tip.  This  derivation 
is,  in  pri  nci  pi  e,  the  same  as  that  given  by  Atkinson  and 
Eshelby  in  consideration  of  the  Broberg  model  ) 

which  deals  with  the  case  of  a  constant  velocity  crack  run¬ 
ning  at  both  ends.  This  provides  the  basis  for  computing 
the  amount  of  energy  released  by  a  running  crack  in  the 
target  material  after  impact. 

The  general  expression  derived  in  Section  2.1  does  not 
depend  on  behavior  of  crack  velocity  with  respect  to  time. 
Subsequent  reduction  is  made  for  uniform  crack  extension 
both  for  a  steady  state  and  non-steady  state  of  stress  ahead 
of  the  crack.  The  general  result  for  both  of  these  cases 
of  uniform  crack  extension  may  be  written  in  terms  of  a  con¬ 
tour  integral  which  Sih  ^  has  shown  to  be  path  independent. 
The  energy  calculation  for  the  problem  of  a  uniformly  ex¬ 
panding  crack  with  a  Dugdale-type  plastic  zone  is  treated 
separately  as  the  general  expression  cannot  be  readily 
applied. 

2.1  GENERAL  DYNAMIC  CASE  FOR  ARBITRARY  CRACK  MOTION 

An  energy  balance  must  be  written  for  a  contour  C  en¬ 
closing  the  tip  of  a  crack  moving  with  a  constant  velocity 
v  as  shown  in  Fi gure  1  . 
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The  fixed  coordinates  are  (X^X^)  and  the  moving  coor¬ 
dinates  attached  to  the  crack  tip  are  (xlSx2)  where 


x i  =  X.  -  vt 6^  ,  (i  =  1,2). 


(1) 


For  a  given  system  of  particles  surrounding  the  crack 
tip  the  increment  of  irreversible  energy  pumped  into  the 
crack  tip  is  given  by  the  relation 

-AaG  =  AE  -  All  (2) 

where  G  is  the  rate  at  which  energy  is  dissipated  by  the 
moving  crack,  AE  is  the  increase  in  total  energy,  and  AU 
is  the  work  done  on  the  system. 

Referring  to  Figure  2,  the  system  of  particles  to  be 
considered  is  included  in  the  area  Aq  bounded  by  C  and  the 
area  A  bounded  by  C  and  C  at  time  t.  At  time  t  +  At  the 
same  particles  are  bounded  by  C  (enclosing  Afl)  and  C  and 
C2  (enclosing  A2).  Note  that  and  C2  are  fixed  contours 
and  C  is  a  moving  contour.  Then 
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Figure  2.  Incremental  Motion  of  Crack  Tip 
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where 


.  u  +  1  Dui  Dui 
5  "  W  2  p  Dt  Dt 


W  +  ^  Puiui 


(5) 


Noting  that  C>  +  C2  coincides  with  C  in  the  limit  and 
dividing  through  by  At,  the  limit  as  At-*-0  is 


-vG  = 


ft  dX.dX2  '  |  5<Xk-t)V6!idSi  - 


Du, 


Dt  '  i 


T,  ds 


or 


v  G  + 


If-  dx  dx 
3 1  1  2 


U(xk,t)vdx2  +  T i u  •  d s ]  .  (6) 


2.2  STEADY  STATE  AND  STATIC  CASE 


For  the  steady  state  case 


i£  =  0 

at 

and 


Du .  3u .  3u . 

_ L  =  _ I  .  v  _ ]. 

Dt  at  axj 


=  -  V 


3u. 

w. 


(7) 


Then 


au. 


Gs.s.  =  t5*(xk)dx2  -  Tn.  33^-  ds] 


(8) 


where 


,  3u.  3u. 

S  (xk)  ■  W  +  2  Py2  75T  • 


(9) 
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For  the  static  case,  the  energy  release  rate  is  also 
given  by  Eq .  (8)  with  the  kinetic  energy  term  being  zero. 
Thus  , 


G 


s  t 


[Wdx2 


3u  . 


Ti  TT  dsl 


(10) 


(12) 

This  result  was  given  by  Rice  .  The  integrals  in  Eqs.  (8) 
and  (10)  are  path  i ndependen t,  a s  shown  by  Sih^  for 
Eq.  (8). 


2.3  BROBERG'S  PROBLEM 


For  the  case  of  a  uniformly  expanding  crack  of  length 
2vt,  H/3t  is  not  zero.  However,  making  use  of  the  dis¬ 
placement  and  stress  solutions  which  are  given  in  Appendix  I 
for  mode  I  opening,  it  may  be  shown  that  in  the  limit  as  the 
area  enclosed  by  the  contour  is  reduced  to  zero  the  following 
expression  holds: 


1  i  m 
OO 


It  dA  =  0  • 

at  o 


Ol) 


Hence, 


vGt  =  lim 
1  OO 


[v£(xk,t)dx 


T  i v  j  d  s  3 


(12) 


Making  use  of  Equations  (5),  (7),  and  (9),  the  preceding 
expression  for  Gj  may  be  written  in  the  form, 


G  T  =  lim 
1  OO 


9u . 


U*(xk)dx2  -  T,  ds]  ♦ 
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1  im  f 

OO  L 


3u  . 

_ l 

at 


2v 


3u  .  3u  . 
1  _ ]_ 

3t  3x 

l 


ds 


-  1  im  \ 
OO  v 


3u 


i 


i  at 


ds 


(13) 


As  was  the  case  for  Equation  (11),  the  second  and  third  integrals 

on  the  right  hand  side  of  Equation  (13)  reduce  to  zero  in  the 

limit  as  OO.  Since  Sih  has  shown  that  the  remaining  inte¬ 
rn 

gral  is  path  independent  ,  the  expression  for  Gj  is  the 
same  as  given  in  Equation  (8)  for  the  steady  state  case.  Thus 

in  terms  of  the  function  F  (s  s  )  which  is  defined  for 

12 

several  cases  of  uniform  crack  extension  in  Appendix  I,  a 
general  expression  which  is  applicable  to  all  cases  of  ur.i 
form  elastic  crack  extension  in  mode  I  that  are  considered 
here  may  be  written 


4 

K+l 


s1(l-sp[4sis!-(l+s“)! 


(14) 


where 


G.  =  lim  Gt 
1  v+0  1 


tt(  k+1  ) 

8y 


k2 


(15) 


with  k  being  the  corresponding  static  stress-intensity 
factor  for  the  particular  problem  under  consideration  and 
with  r 


< 


3-4v  Plane  Strain 

(3-v) (1 +v)  >  P1 ane  Stress. 


(16) 
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2.4  PLASTIC  ENERGY  DISSIPATION 


Consideration  of  the  effects  of  ductility  upon  dynamic 

(1  3) 

crack  extension  may  be  made  by  use  of  the  Dugdale 
hypothesis  which  is  explained  in  detail  in  Section  IV.  For 
steady  state  problems  where  the  plastic  zone  size  remained 
fixed  with  time,  the  energy  release  rate  (or  plastic  energy 
dissipation)  may  be  calculated  by  direct  application  of 
Equation  (8)  to  an  equivalent  elastic  problem. 

Consider,  however,  the  motion  of  an  expanding  crack  of 
length  2c  where  c  is  a  function  of  time  and  with  a  thin 

Dugdale-type  plastic  zone  located  ahead  of  the  crack  tip. 

(14) 

As  shown  by  Atkinson  (see  Section  V)  for  uniform  ex¬ 

tension,  this  plastic  zone  must  be  expanding  with  time  in 
order  that  the  Dugdale  hypothesis  of  finite  stress  at  the 
end  of  the  plastic  zone  be  satisfied.  Therefore,  the  re¬ 
sults  of  the  first  part  of  this  chapter  cannot  be  used  here 
even  for  the  uniformly  expanding  crack  since  the  contour  C, 
which  must  enclose  the  entire  region  where  energy  dissipa¬ 
tion  takes  place,  was  assumed  to  be  fixed  in  shape  and  size. 
Obviously,  an  expanding  plastic  zone  size  would  demand  an 
expanding  contour  and  therefore  a  separate  analysis  for 
this  case  is  made  here. 

Figure  3  shows  the  conf i guation  of  the  crack  surface 
and  plastic  zone  boundary  at  times  t  and  t+At  during  which 
time  the  crack  edge  extends  by  an  amount  Ac  =  v(t)At  and 
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Figure  4. 


the  plastic  zone  extends  by  an  amount  Aa  =  B ( t )  A t  (Aaj>Ac). 

The  expression  for  the  increment  of  energy  released  on  the 
right-hand  side  of  the  crack  is  then  equal  to  the  plastic 
energy  dissipated  in  the  movement  shown. 

R. calling  the  issumption  of  a  thin  plastic  strip,  the 
elastic-plastic  problem  being  considered  may  be  converted 
to  the  elastic  problem  of  a  crack  of  length  2a  with  tensile 
stress  equal  to  the  yield  stress  Y  of  the  material  acting 
on  the  region  c<x<a  of  the  crack  surface  (see  Figure  4). 

The  negative  work  done  by  these  tractions  is  then  equal  to 
the  plastic  energy  dissipation  and,  for  the  case  of  per¬ 
fectly  plastic  material,  where  Y  is  constant,  the  resulting 
expression  for  energy  released  during  the  interval  f  time 
At  i  s 

X."af(t)  3U, 

G j AC  =  -AUp  =  2Y  At  d^  (17) 

Xj=c(t) 

where  Gj  is  the  energy  release  rate  per  unit  crack  extension, 

AU  is  the  work  done  on  the  elastic  plastic  boundary,  and 
P 

u2  is  the  displacement  of  the  crack  surface  in  the  x^direc- 
ti  on . 

Now  introduce  the  dimensionless  time  variable  t  as 

t  =  Cjt  (18) 

and  normalize  the  crack  velocity  v  and  plastic  zone  velocity 
6  with  respect  to  c ^  so  that  the  following  normalized  crack 
velocities  are  defined: 

14 


I 


(19) 


Then,  the  result  in  the  limit  as  At-*-0  is 


xra,(t)  9u 

Gib  =  2V  j  37idxi 

X  =c  (t  ) 

1 


(20) 


The  calculation  of  Gj  for  uniform  crack  extension,  u.with 
b  and  constant,  is  carried  out  in  Section  V. 


SECTION  III 


SUDDEN  APPEARANCE  OF  A  CRACK 

The  response  of  a  solid  medium  to  the  sudden  applica¬ 
tion  of  loads  or  the  sudden  appearance  of  geometric  defects 
has  engendered  considerable  interest  in  the  past.  However, 
it  has  been  only  recently  that  any  consideration 
has  been  given  to  the  propagation  of  transient  elastic  waves 
in  bodies  containing  mechanical  flaws  or  cracks.  Due  to  the 
crack  geometry  this  class  of  problems  involves  mixed  edge 
conditions  where  both  stress  and  displacement  must  be  pre¬ 
scribed  on  a  particular  boundary.  Such  conditions  can  give 
rise  to  singular  stresses. 

An  important  case  associated  with  the  problem  of  pro¬ 
jectile  penetration  is  that  of  a  crack  appearing  suddenly 
in  a  stressed  medium.  The  associated  disturbance  due  to  the 
■  c  may  be  studied  in  two  separate  regions:  the  local  re¬ 
gion  near  the  crack  edge  (near-field)  or  the  region  close 
to  the  wave  front  propagating  away  from  the  crack  (far- 
field).  In  fracture  mechanics,  much  of  the  interest  is 
centered  on  the  near-field  solution  which  describes  the 
condition  of  possible  crack  extension.  The  far-field  solu¬ 
tion  may  be  of  use  in  determining  the  effects  of  the  crack 
on  other  parts  of  the  structure.  The  necessity  of  descri¬ 
bing  the  two  solutions  separately  is  due  to  the  adopted 
mathematical  technique  for  solving  crack  problems  which 
generally  yields  a  tractable  solution  in  only  one  of  these 
regions.  ,, 


A  number  of  crack  problems  of  this  type  have  been 
solved  for  the  case  of  anti-plane  shear  (see,  for  example, 
Kostrov  '  '  );  however,  the  mathematical  techniques  employed 

are  so  far  inapplicable  to  problems  of  in-plane  loading. 

For  in-plane  loading,  the  only  crack  problem  that  has  been 
solved  so  far  is  concerned  with  the  case  of  a  semi - i nf i n i te 
crack  suddenly  appearing  in  a  stretched  elastic  plate. 

Maue  ('^and  Baker  have  considered  this  problem  using 

the  Weiner-Hopf  technique.  The  solution  obtained  is  un¬ 
bounded  as  t-*°°  and  does  not  permit  meaningful  comparison 
with  a  static  case;  however,  the  near-field  stress  solution 
is  valid  near  the  tip  of  a  finite  crack  for  the  initial 
period  of  time  before  interaction  with  stress  waves  emana¬ 
ting  from  the  other  crack  tip  occurs. 

In  this  section,  the  complete  near-field  solution  to 
the  problem  of  a  finite  line  crack  subjected  to  transient 
in-plane  loading  will  be  considered.  The  transient  response 
of  the  material  in  the  neighborhood  of  the  crack  tip  will  be 
obtained  for  the  case  of  sudden  application  of  self-equili¬ 
brating  tractions  to  the  crack  surface.  The  crack  is  con¬ 
tained  in  an  initially  unstressed  infinite  elastic  plate 
and  is  oriented  as  shown  in  Figure  5.  For  mathematical 
simplicity,  the  load  is  separable  into  a  spatial  and  temporal 
function  where  the  spatial  function  is  even  in  x  .  Symmetry 
conditions  then  allow  the  problem  to  be  solved  in  the  upper 
plane  only  .  The  special  case  where  sudden  uniform 
pressure,  aH(t),  is  applied  is  treated  in  detail.  Here, 
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Self -Equilibrating  Tractions 
q(x,)  f(t) 


Figure  5. 


Crack  Geometry  and  Stress  Components 


H(t)  is  the  Heaviside  step  function.  The  addition  of  this  pro 
blem  to  the  case  of  an  infinite  uncracked  plate  under  uni¬ 
form  tensile  stress  a  corresponds  to  the  problem  of  the 
sudden  appearance  of  a  crack  in  a  stressed  plate. 

The  analysis  is  done  in  two  major  steps.  The  problem 
is  first  formulated  with  the  aid  of  Fourier  and  Laplace 
transforms  and  reduced  to  a  Fredholm  integral  equation  in 
terms  of  the  Laplace  transform  variable.  The  Laplace  trans¬ 
forms  of  the  stress  components  are  then  inverted  by  a  com¬ 
bination  of  numerical  means  and  application  of  the  Cagniard 
inversion  technique. 

3.1  FORMULATION  OF  PROBLEM 


For  plane  elastodynamic  problems,  the  displacement  com¬ 
ponents  may  be  expressed  in  terms  of  two  scalar  potentials, 
^>(xi  *X2  *t)  anc*  4'(x1>*2>t)  so  that 


94)  9i|) 

9  x  9x 

2  i 


(21) 


where  4>  and  ip  satisfy  the  second-order  partial  differential 
equations 


V24> 


1_  9^ 

C2  9t 2 
l 


I_  9 

2  , 
c  9t 2 
2 


(22) 


with  V2  being  the  Laplacian  operator.  The  dilational  and 
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shear  wave  velocities  are,  respecti vely  , 


Cj  =  [(X+2u)/p]*  ,  c2  =  ( m/ p ) 2 


(23) 


where  \  and  p  are  the  Lamd  constants. 


Applying  the  strain  displacement  and  constitutive  equa¬ 
tions  for  a  homogeneous  isotropic  elastic  media, the  stresses 
may  be  written  in  terms  of  d>  and  \ p  as, 


o 


i  i 


-pc*V2<f>  +  2pc* 


+ 

3x2  hXl 


a 

2  2 


-pc  j  V24> 


3jJj_  dtp  ] 

3x2  '  dxj 


1  2 


d2ip  _  3^1 
k3x2  3x2 


+ 


o  32  $ 
3x  j  3x 


(24) 


If  the  material  is  initially  unstressed,  the  stress 
and  displacement  field  due  to  sudden  application  of  normal 
tractions  to  the  crack  surface  may  be  found  by  solving  the 
preceding  field  equations  subject  to  zero  initial  conditions 
and  the  following  boundary  conditions  at  x2  =  0: 


a2  2 ( x i ’ 0 , t ) 

ai2 (XJ ,°,t) 
u2(Xj  ,0,t) 


In  addi  ti  on,  the 


=  -oq ( x  j ) f ( t )  ,  for  i x  j  j  < a 

=  0  ,  for  0<|x1  |<°°  (25 ) 

=  0  ,  for  | x  | >a  . 

condition  on  displacement  at  infinity  is 
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I 

lim  [ui(xi,x2,t)  ,  u2(xi»x2,t)]  =  0  . 

X  2  +  X  '  +  00 

1  2 

The  parameter  a  is  a  constant  with  the  dimension  of  stress 
and  q(x^)  is  restricted  to  functions  that  are  even  in  xt; 
that  is, 

q(Xj)  =  q(-xx )  . 

3.2  DUAL  INTEGRAL  EQUATIONS 

Recalling  that  the  initial  conditions  are  zero,  the 
Laplace  transform  may  be  applied  to  Equation  (22)  with  the  re¬ 
sult 

V24>*  =  (p/Cj)2^*  ,  vV  =  (p/c2)2^*  (26) 

where  the  Laplace  transform  pair  is  defined  by  the  equations 


00 

f*(p)  =  |  f (t)exp(-pt)dt 
0 

f(t)  =  2~j  |  f*(p)exp(pt)dp 
where  the  second  integral  is  over  the  Bromwich  path 


In  order  to  reduce  Equation  (26)  to  ordinary  differ¬ 
ential  equations,  the  Fourier  cosine  and  sine  transforms  will 
be  applied.  The  Fourier  cosine  transform  pair  is 


f  (a) 


00 

* 

f(x  )cos(ax  )dx 

i  1  1 


f(xt)  =  f  |  f(o)cos(axj )da 
0 
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and  the  Fourier  sine  transform 


f(ct)  = 


f (x, )sin(ax  )dxx  ,  f(xj)  =  -  f (a)sin(ax1  )da 


The  application  of  these  transforms  will  depend  on  whether 
the  function  under  consideration  is  even  or  odd  in  x  . 

i 

By  consideration  of  the  symmetry  properties  of  the 
boundary  condi tions, i t  can  be  shown  that  for  the  case  con¬ 
sidered  here, 


<f>(xlSx2,t)  =  cH-Xj  ,x2  ,t)  ,  iMxlSx2,t)  =  -4»(-x1,x2t) 

t'i 

and  the  solution  may  be  considered  in  the  first  quadrant  of 
the  plane. 

Then,  the  Fourier  cosine  transform  is  applied  to  the 
function  <p*  and  the  Fourier  sine  transform  is  applied  to 
ip*  in  Equation  (26),  the  result  being, 


-  IV  +  (p/c  yjr  -  0 

dx2  1 

^  -  [a2  +  (p/c2)2]J*  =  0  . 

dx22 


Solutions  to  Equation  (27)  that  satisfy  the  regularity 
conditions  at  infinity  are 


<P*  ■  |  A(a,p)exp(-y1x2) 

A 

4>*  =  \  B(a,p)exp(-Y2x2 ) 


(28) 
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In  order  to  put  the  dual  integral  equations  in  a  form 
amenable  to  solution,  the  first  of  Equation  (31)  is  integrated 


with  respect  to  x:  over  the  interval  (O^). 


Tne  result  is 


7  ®Q(x  )f*(p) 

g(ak)D(a,p)si  n(ax ,)da  =  - 

1  2 

0  pp 

00 

[  D(a,p)cos  (ax  i  )da  =  0  |xj>a 
0 


x  3  |  c  a 


(32) 


where 


Q(xx) 


■  i 

q(s)ds  . 

0 


Following  the  usual  argument  for  solution  of  dual  inte¬ 
gral  equations  ,  a  displacement  related  function  h*(x  ,p) 
is  defined  by  the  equation. 


h*(xl,p) 


00 

* 

D(a,p)cos (ax 


)  da  . 


(33) 


Then,  noting  that  the  second  of  Equation  (32)  requires  that 
h*(x  ,p)  be  zero  for  lx  | > a ,  the  Fourier  inversion  theorem 

i  i 

leads  to  the  expression 


D(a,p) 


a 

h* ( x  , p) cos (ax i )da  . 

0 


(34) 


The  function  h*(x  ,p)  is  to  be  constructed  to  possess  the 
proper  asymptotic  behavior  at  the  crack  border,  so  that  the 
near-field  displacement  is  proportional  to  the  square  root 
of  the  distance  from  the  crack  edge.  To  this  end,  let 
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I 


h 

\ 


h* ( x  t  ,p ) 


U  (x  , p)xdx 

-/x  2  -  x  * 


0< x j <a  . 


(35) 


Substitution  of  Equation  (35)  into  Equation  (34)  and  application  of 
the  identity 


i 

f 


cos ( ax 


0 


/i 2  - x  2 


(ax  ) 


(36) 


where  J  denotes  the  zero  order  Bessel  function  of  the  first 
o 


kind,  leads  to  the  expression 


D(a,p)  =  | 


a 

r 


U  (x  ,p)  J  Q (ax )xdx  . 


(37) 


In  order  to  rewrite  the  first  of  Equation  (32)  as  an 
Abel's  integral  eouation  for  U*(x,p)  a  function  W(ak)  is 
introduced  such  that, 


W(ak)  =  (l-cj)-g(ak) 


(38) 


By  virtue  of  Equations  (37)  and  (38),  the  governing  integral 
equation  takes  the  form 

ar  *  1  ?  °Q(x,)f*(p) 

(1-c2)  [U  ( x  ,  p )  x  J0(ax)sin(axl)da]dx  =  -  -  + 


J 

0 


0 


PP 


U  (  T  ,  p  )  X 


W(ak  )  J  0  ( ax  )  s  i  n  (ax  j  )da]dx  |xj<a  . 

(39) 
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L 


Using  the  identity. 


CKx  j  <T 


J  Q  (  ax  )  s  i  n  (  ax  ^ )  da  =  <{ 


(x*  -  t2  )"^  x  > T 


the  desired  form  of  Abel's  integral  equation  is  obtained 


(1-cJ) 


o 


u 


pp‘ 


a  00 

+  j  [  U  *  (  T  ,  p  )  T 
0 


W(ak)J0(ai)sin(ax1)da]dT 


x  <a 

1  1 


Inversion  of  the  preceding  equation  yields 


(1-cpu‘O.p)  =!{-'— 

0  /rTT*T  l  pp 


+  j  [u*(c,p)c 
0 


W(ak)JQ (ac)cos(axi )adct]dc |  dxt 


t  <  a 


(40) 


Now,  introduce  the  non-dimensional  variables 


K  =  i/a  ,  n  •-  C/a  ,  s  =  aa  ,  z  =  x^a 


and  define 
<  =  k/a 


S ( z )  =  W(za) 


A*(e,K)  -  £  n-cM  pP  1/7  u  ^“-p’  • 

2  0  of*(p) 
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Then,  making  use  of  Equation  (36),  Equation  (40)  may  be  written  in 
the  form  of  a  Fredholm  integral  equation  of  the  second  kind. 


A  U,k)- 


A  (n.<)KU,n)dn  =  - 


0 


B(z)dz 

/C2-z2 


(41  ) 


whose  kernal ,  being  symmetric  in  £  and  n,  is 


KU,n)  = 


/fn 


n-c;i  o 


sW(sk)  J0  (sn)Jo  (s£)ds  ,  0<f;£l  ,  0<n<l  . 


For  rapid  convergence  of  the  infinite  integral,  define 
a  function  d(z)  as 

d(z)  =  W(z)  + 

z  2  +  E2 

where,  in  order  that  d(z)  be  of  order  (z)-6  for  large  z, 

H  and  E2  are  chosen  as 

H  =  1  ( 3 c g  -  4cJ  +  3) 

r,  [5c5  -  6c;  +  2cJ  +  1] 

E  =  8H 

Since 

oo 

1  7  TI  ra  Jo(sn)J0(sc)sds  =  !L  i„  (^)K„(^).  0<em  , 

5  (  ClK  )  +E  K 

the  expression  for  the  kernal  in  the  Fredholm  integral 
equation  becomes 
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K  ( € ,  n )  = 


/Cn 


(T-Cfl)  L 


^  lt(P)K„(^)  ♦ 


0  '  < 


sd(sK)J0 (sn)J0 (sC)ds 
0<C<n 


where  Ifl  and  Kfl  are  the  zero-order  modified  Bessel's  func¬ 
tions  of  the  first  and  second  kind  respectively. 


Making  use  of  Equation  (37)  and  the  definition  of  A*  the 
unknown  function  D(a,p)  may  be  written  in  terms  of  the 
solution  to  the  Fredholm  integral  equation  as, 


D(ct,p) 


aa2  f * ( p ) 

1 

' 

a*U,k)1 

(l-cj)pp2  j 

) 

L  /c 

J  0 (aa£ )£d£  . 


This  expression  may  be  integrated  by  parts  with  the  result 


D(a,p)  =  -aaf*M  ...  (a*(i  ,<)j  (aa) 
(l-cMpp2cx  V 


d 

"a*U,k)‘ 

d£ 

0 

L/c  J 

J  (aaC)CdC 


(42) 


from  which  the  stresses  and  displacements  may  be  determined, 


3.3  TIME  DEPENDENT  STRESS  FIELD 

The  integral  expressions  for  the  Laplace  transforms  of 
the  potentials  $  and  ip  are  now  completely  determined  and, 
making  use  of  Equation  (24),  corresponding  expressions  for  the 
Laplace  transforms  of  the  dynamic  stresses  may  be  obtained. 
It  remains  then  to  take  the  inverse  Laplace  transforms. 

This  may  be  accomplished  by  applying  the  Cagni ard-DeHoop 
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inversion  technique  ^  ^ .  Due  to  the  complexity  of  cal¬ 
culation,  the  detailed  procedure  will  be  given  only  for  the 

first  stress  invariant,  a  +  o 

11  2  2 

From  Equations  (24)  and  (26)  the  Laplace  transform  of 

o  +o  is 
11  2  2 

o*i  +  o*2  =  -2pp2  ( 1  -c]  )<j>*  (43  ) 


The  part  of  the  solution  of  interest  here  is  the  singu¬ 
lar  solution  near  the  crack  tip.  For  this  portion  of  the 
solution  it  is  sufficient  that  only  the  first  term  on  the 
right  hand  side  of  Equation  (42)  be  retained.  The  integral 
expression  in  D(a,p)  is  associated  with  terms  that  remain 
finite  at  the  crack  tip. 

Putting  <p*  in  Equation  (30)  into  Equation  (43),  the  near-field 

solution  for  o,,  +  a,,  in  the  transformed  plane  is 
1112 


0  ,  +0 
11  2  2 


=  -2oaf*(p)A*(l  ,tc) 


ay 


’  [a2  +  1  (p/c2)2] 


•  J 1 (aa ) exp ( -y ^  ) cos (ax  i  )da 


(44) 


Changing  the  variable  of  integration  to  a=  wp  and  noting 
that  J  (awp)  is  an  odd  function  with  respect  to  w.  Equation  (44) 
becomes 


a 


★ 

1  1 


+ 


oo 

-craf*(p)A*(  1  ,<)p 

-a 


T(w)Jj(apw)exp[-p(  /w2+c  j5'x2 
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in  which  T(w)  stands  for 


T(w)  = 


w2  +  ( 1 / 2c2  ) 


w/w2+( 1 / C  2  ) 


Making  use  of  the  identity 


(21  ) 


J,(x)  =  i 

1  '  '  7T 


exp(-ixcosto)cos(co)dcu 


Equation  (44)  further  simplifies  to 

-  f*(p)A*(l  ,k) 


*  * 

O  +0 
11  2  2 


cos  (to)du)j 


T  ( w ) 


*  exp [-p(/w2+c‘2x2-iwX)] 
The  variable  X  is  given  as 
X  =  x j-  acosui  . 


dw 


(45) 


To  evaluate  Equation  (45),  consider  first  the  function  I 
defined  by  the  relation 


I 


7T 

r 


costodw 


T(w)exp[-p(  /w2+c  ~’5  x2  -  i  wX )  ]dw 


T4fi ) 


This  function  is  in  a  form  suitable  for  evaluation  through 
the  Cagni ard-DeHoop  method.  The  objective  then  is  to  trans¬ 
form  the  integral  on  the  right-hand  side  of  Equation  (46)  into 
a  recognizable  Laplace  transform.  This  may  be  done  by  making 
the  change  in  variables, 


31 


I  =  v  w 2  +  c-2  x,  -  iwX  . 

1  2 

Solving  for  w  as  a  function  of  t  yields 

X 

w  =  +/?2"-  (R/c,  )2  —  +  —  (47) 

1  R2  R2 

in  which  R2  =  X2  +  y2. 


The  remainder  of  the  procedure  consists  of  treating  w 
as  a  complex  variable  and  locating  a  path  in  the  w-plane 
along  which  t  is  a  positive  number.  Referring  to  Figure  6 
and  noting  that  the  integrand  in  Equation  (46)  possesses  branch 
points  at  + i c  ,  cuts  are  made  in  the  w-plane  as  shown  in 
order  that  the  f unc ti on Vw2  +  ( 1  /c  1  )2  be  single  valued  and 
have  a  real  part  that  is  positive  everywhere  in  the  w-plane. 
If  t  is  restricted  to  real  non-negative  values.  Equation  (47) 
is  the  equation  of  a  hyperbola  in  the  w-plane.  If  X>0  (X<0) 
the  hyperbola  will  be  in  the  upper  (lower)  half  of  the  w- 
plane.  The  path  of  integration  in  Equation  (46)  can  be  converted 
from  the  real  axis  to  the  aforementioned  hyperbola  by  making 
use  of  the  Cauchy  integral  theorem  and  Jordan's  Lemma. 

Thus,  after  changing  the  variable  of  integration  to  t  by 
means  of  Equation  (47),  it  is  found  that 


coswdoj 


r  -■■  /  +\  dw  + 
[T(w  )  37- 


-  T(w)  37— ]exp(-px)dT 


k/c 


(48) 


where  the  (+)  and  (-)  superscripts  refer  to  the  right-  (left-) 
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Figure  6,  Path  of  Integration  in  w-plane. 


hand  legs  of  the  hyperbola.  Noting  that  the  inner  integral 
is  in  the  form  of  a  Laplace  transform  it  is  a  matter  of 
inspection  to  write  the  inverse  Laplace  transform  of  I*  as. 


i  (t) 


{  H(t-  §-)  CT(w+) 

0  1 


T(w‘)  ]coso)daj  (49) 


where  H(t)  is  the  Heaviside  step  function.  Referring  to 
Equations  (45)  and  (46)  and  making  use  of  the  convolution 
theorem  for  the  Laplace  transform,  the  stress  invariant  is 
obtained  to  be 


a  +o  =  -  — —  m(t-T)I(x)dT 

1  1  2  2  7T  '  ' 


where 


m(t)  =  L" * [pf* ( p  )A* ( 1  ,k  )  ] 


with  L  1  being  the  inverse  Laplace  transform  operator. 


In  order  to  determine  a  +  a,  in  the  region  near  the 
crack  tip  it  is  necessary  to  extract  the  singular  portion 
of  I  ( t )  from  Equation  (49).  Without  going  into  mathematical 
details,  it  is  easy  to  show  that,  as  long  as  c2t  is  suffi¬ 
ciently  greater  than  r,  the  radial  distance  measured  from 
the  crack  tip,  I  is  independent  of  time,  i.e., 


(q-z)/a2-q2 


(q-z)/a2-q‘ 


L 


with 


2  =  x>  +  ix2 


2  =  X  ,  -  1  X 

1  2 


and  where  the  transformation  of  variable, 
q  =  acosu 

has  been  made. 

Making  use  of  Chapter  4  of  Mus khel i shvi 1 1 i  ^ ^ ,  the 
evaluation  of  the  above  integrals  near  the  crack  tip  leads 
to  the  following  expression  for  the  singular  portion  of  I, 


I  = 


2m' 
/2T r 


cos(e/2)  +  • 


where  r  and  6  are  polar  coordinates  with  the  origin  attached 
to  the  crack  tip  as  shown  in  Figure  7. 

Then,  referring  to  Equation  50,  it  follows  that  the  near 
field  solution  for  the  stress  invariant  is 


a  +  a, , 
11  2  2 


2a/a 

/2T 


M(t)cos(6/2)  + 


(52) 


where 


M(t)  = 


m(t)dt  . 


(53) 


In  a  manner  analagous  to  the  static  theory  of  brittle 
fracture,  a  dynamic  stress  intensity  factor,  k^t)  may  be 
defined  ^ 


Si 
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Outgoing  Waves 


\ 

s 


Figure 


7.  Wave  Fronts  Emanating  from  Crack  Tip 


k, (t)  =  M(t)  a  . 


(54) 


With  this  definition  the  individual  stress  components  may 
be  evaluated,  the  result  being 

Mt) 

o  =  -  cos(e/2)[l-sin(6/2)sin(30/2)]  +  ••• 

11  v7r 

k,  (t) 

a,,  =  - cos (0/ 2) [1 +s i n (6/2 ) si n ( 36/ 2 ) ]  +  ••• 

22  /Zr 

k  j  ( t ) 

a  =  -  cos(0/2)sin(0/2)cos(30/2)  +  •••  (55) 

v7r 

Thus,  the  stress  field  near  the  crack  tip  has  the  same 
spatial  distribution  for  the  dynamic  case  as  for  the  static 
case  the  only  difference  being  that  the  intensity  of  the 
field  is  a  function  of  time.  The  time  dependent  stress- 
intensity  factor  k:(t)  is  related  to  the  function  m(t)  which, 
noting  Equation  (51)  depends  ultimately  on  the  function  A  (1,k). 
Thus,  in  order  to  obtain  the  stress- i ntens i ty  factor  for  a 
given  time  dependent  load  applied  to  the  crack  surface,  it 
is  necessary  to  first  numerically  solve  the  integral  Equa¬ 
tion  (41)  for  A* ( 1 , k )  for  various  values  of  the  Laplace 
transform  variable  p.  Inversion  of  the  Laplace  transform 
numerically  then  enables  the  evaluation  of  M(t)  and  thus 
kj(t).  This  procedure  is  described  in  Appendix  II. 

Numerical  results  have  been  obtained  for  the  case  of  a 
uniform  load  a  suddenly  applied  to  the  crack  surface.  That 
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is,  jq ( x  i  )  is  a  constant,  and  f(t)  is  the  Heaviside  step 
function  H(t).  By  superposi tion,  this  type  of  loading  cor¬ 
responds  to  the  sudden  appearance  of  a  crack  in  a  stressed 
plate.  The  solution  to  the  Fredholm  integral  equation, 
A*(1,k)  for  this  case  is  plotted  in  Figure  8.  In  Figure 
9  the  dynamic  stress-intensity  factor,  k^t),  normalized 
with  respect  to  the  corresponding  static  value  a/a ,  is 
plotted  as  a  function  of  c2t/a.  The  ratio  of  shear  wave 
speed,  c2,  to  dilatational  speed  c,  that  was  used  is  0.542 
the  value  for  steel.  The  significant  feature  of  the  solu¬ 
tion  is  that  the  stress- intensity  factor  reaches  a  peak 
greater  than  the  static  value  and  subsequently  oscillates 
about  that  value  with  decreasing  amplitude  until,  in  the 
limit  as  t-*°°,  k}  -  c/a.  When  k,  reaches  its  maximum  value, 
c2t/a  is  approximately  equal  to  3.0.  For  a  crack  where  a 
is  equal  to  one  inch,  this  corresponds  to  an  elapsed  time  of 
2.4  x  10  5  seconds. 

The  method  used  here  to  determine  the  asymptotic  be¬ 
havior  of  the  near- 'field  stress  solution  requires  only  that 
the  stress  wave  boundary  be  well  beyond  the  region  of  in¬ 
terest.  Thus,  if  the  stress  field  is  to  be  valid  within  a 
distance  c  of  the  crack  tip,  the  condition  on  the  time 
elapsed  is  c2t>>t:.  As  long  as  c  is  kept  small  enough,  the 
solution  is  valid  for  all  time  (except  very  near  t  -  0). 
This  condition  ^n  G  is,  of  course,  in  addition  to  the  usual 
static  restriction  that  e<<a. 
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Figure  8.  Solution  to  the  Fredholm  Integral  Equation. 


SECTION  IV 


PLASTICITY  MODELS  OF  STEADY  STATE  CRACK  EXTENSION 

It  is  well  known  that  upon  reaching  a  certain  critical 
stress  or  strain  level  the  material  in  a  region  local  to  the 
crack  tip  will  separate,  caus i ng  crack  propagation  in  a 
manner  determined  by  the  character  of  the  stress  field  caused 
by  impact  and  the  structure  of  the  target  material.  In  ma¬ 
terials  without  a  substantial  amount  of  yielding,  the  crack 
will  normally  start  its  motion  abruptly  and  then  very  ouickly 
settle  down  to  a  uniform  speed.  To  estimate  the  amount  of 
energy  carried  by  a  running  crack,  it  is  essential  to  have 
a  mathematical  description  of  the  stress  state  near  the  crack. 
This  is  important  in  determining  the  condition  required  for 
cessation  of  crack  motion. 

For  static  problems,  plasticity  effects  have  been  in¬ 
troduced  into  fracture  analysis  with  a  satisfactory  amount 
of  success  (see  a  review  article  by  Rice  for  example). 

However,  very  little  analysis  has  been  done  which  deals 
with  plastic  behavior  of  material  in  the  neighborhood  of 
propagating  cracks.  For  the  few  elastic  solutions  to  dynamic 
crack  problems  which  exist,  the  introduction  of  plastic 
effects  causes  complexities  which,  with  the  exception  of 
the  Dugdale  model,  have  been  beyond  resolution. 

The  elastic  crack  propagation  problems  that  have  been 
studied  are  in  two  basic  configurations.  The  solutions  of 
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/ ^  r  \  (  0  £  \ 

Yoft'6  c  and  Craggs  '  '  are  representative  of  the  first 

type  (steady  state).  Yoff6  considered  the  problem  of  a 
crack  of  fixed  length  moving  through  an  infinite  elastic 
continuum  with  a  constant  velocity.  Craggs  considered  a 
similar  proDlem  where  a  semi -i nf i ni te  crack  propagates  at 
a  constant  velocity  and  upon  which  tractions  are  applied  to 
a  segment  of  the  surface  which  moves  with  the  crack  tip. 

Broberg's  analysis  of  the  physically  more  realistic 

configuration  of  a  crack  whose  tips  are  expanding  away  from 
the  origin  at  a  constant  rate  is  r epres enta t i ve  of  the 
second  type.  The  solution  obtained  by  Broberg  was  also  ob¬ 
tained  by  Craggs  who  made  use  of  the  dynamic  similarity 

possessed  by  solution. 

4.1  THE  DUGDALE  MODEL 

One  available  tool  for  extension  of  these  configurations 
to  include  plastic  behavior  is  the  Dugdale  model.  The 

Dugdale  hypothesis,  proposed  for  the  case  of  plane  stress 
yielding  ahead  of  the  crack  tips,  postulates  that  the  yielded 
zone  is  a  narrow  region  in  the  plane  of  the  crack  (see  Fig¬ 
ure  10).  For  a  stationary  crack  of  length  2c,  the  length 
of  this  zone  is  calculated  by  solving  the  elastic  problems 
of  a  crack  of  length  2a  with  remotely  applied  tensile  stress 
P  and  with  tensile  stress  Y  distributed  over  a  length  (a-c) 
of  the  crack  which  represents  the  boundary  of  the  plastic 
zone  (see  Figure  11).  A  relationship  between  the  geometry 
of  the  crack  and  plastic  zone,  the  applied  load, and  the 
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yield  stress  Y  is  obtained  through  the  requirement  that 
stress  be  finite  at  the  end  of  the  plastic  zone.  The  ex¬ 
pression  thus  obtained  is 

I  “  cos(j  y)  .  (56) 

This  result  was  compared  favorably  with  experimental  measure¬ 
ments  by  Dugdale. 

Subsequently,  the  Dugdale  model  has  been  incorporated 

into  the  case  of  a  propagating  semi-infinite  crack  by  Goodier 

and  Field  the  Yoff£  model  by  Kanninen  et  al  and 

(14  1 

the  Broberg  model  by  Atkinson  v  .  In  this  section  and 
Section  V,  these  models  will  be  examined  with  special  atten¬ 
tion  directed  toward  the  energy  dissipation  associated  with 
each . 

4.2  STEADY  STATE  CRACK  PROPAGATION  PROBLEMS 

For  steady  state  problems,  with  a  disturbance  moving  to 

the  right  with  velocity  v  in  an  infinite  or  s emi - i nf i ni te 

medium, the  Galilean  transformation  of  coordinates,  x-  =  X  - 

1  i 

vt5^,  allows  the  reduction  of  the  number  of  independent 
variables  from  three  to  two.  Thus,  many  of  the  powerful 
tools  of  two-dimensional  analysis  may  be  brought  to  bear 
on  the  problem.  In  parti  cul  ar,  the  comp'lex  variable  technique 
of  Muskhel ishvill i  is  of  great  use  in  dealing  with  steady 

(in)  (ill 

state  crack  propagation  problems  (see  References  v  ,  ,  and 

(321 

v  ' ) .  Thus,  for  a  crack  moving  with  velocity  v,  the  stress 
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components  of  dynamic  plane  elasticity  may  be  expressed  in 

terms  of  two  functions  (^1)  ,  <*>.(  z.),  (j=l,2),  of  the  complex 

3  3 

variables  z  .  as 
3 

aM  =  -2Re[Sj  +  1(  1  -s'  )  t ,  (z , )  +  [( 1+s’  )*2  ( 2;  )  ]  (57) 

°22  =  O+spReC^^Zj)  +  <i>2(z2)]  (58) 

n+sM 

°i2  =  2Im[si$l(zi)  +  -  - . -  *2(z2)]  (59) 

2 

where  sj  and  s^  are  defined  in  Appendix  I  and, 

Zj  =  X1  +  iSjX2  ,  j  =  1 ,2  (60) 

with  (x1»x2)  being  the  coordinates  in  a  frame  of  reference 
attached  to  the  crack. 

Similarly,  the  displacement  field  is 

u,  =  -  J-  ReO^Zj)  +  l(l+s^)<}>2(z2)]  (61) 

u2  =  l  Im[s  (z,  )  +  MZ2)]  (62) 

in  which 

VZj>  =  1  *j(zj,dzj  •  j  5  ''2  •  (63) 

Due  to  simplicity  of  analysis  then,  most  crack  propa¬ 
gation  problems  that  have  been  solved  are  steady  state.  The 
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two  examples  considered  in  this  section  are  a  semi - i nf i ni te 
crack  with  self-equilibrating  concentrated  loads  Q  applied 
at  a  fixed  distance  c  from  the  crack  tip  and  a  finite  crack 
of  fixed  length  2c.  Both  cracks  propagate  with  constant 
velocity  v. 

4.3  SEMI-INFINITE  CRACK 


Considered  here  is  the  case  of  a  semi -infinite  crack 
whose  motion  is  maintained  by  concentrated  loads,  Q,  at  a 
distance  c  from  the  crack  tip.  The  geometry  is  shown  for 
the  elastic  case  in  Figure  12  and  for  the  plastic  case 
in  Figure  13. 


The  stresses  and  displacements  for  the  elastic  model 
may  be  obtained  from  Appendix  I.  Makinu  use  of  Equation  ( 14) , 
the  elastic  energy  release  rate  for  plane  stress  is 


Gj  (l+v)s  j (1-s*) 

* 1  =  4si$2  -  (1+s*)2 


(64) 


★ 

where  Gj  is  the  corresponding  static  value  of  the  energy 
release  rate  and  s-  ( j  =  1  , 2 )  is  defined  by  Equation  (I-f).  For 

J 

★ 

the  semi -infinite  model  Gj  is 


★ 


1  i  m  G 
v+0 


(65) 


★ 

The  ratio  G j / G ^  is  plotted  in  Figure  18.  This  result  is 

also  valid  for  the  case  of  the  fixed  length  finite  crack 

★ 

(with  appropriate  change  in  Gj). 
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Figure  12.  Semi -infinite  Crack  Loaded  with  Concentrated 
Forces . 


Q 


*2 


Plastic  Zone  (cr22= Y) 


Using  the  Dugdale  hypothesis  a  thin  plastic  zone  ahead 
of  the  crack  may  be  represented  by  an  extension  of  the  crack 
by  a  length  equal  to  that  of  the  plastic  zone.  Over  this 
added  increment  of  crack  length  tensile  stress,  a22  =  V 
equal  to  the  yield  strength  of  the  material  is  applied.  This 
is  shown  in  Figure  14. 

The  incorporation  of  plasticity  into  the  semi -i nf i ni te 
model  may  then  be  accomplished  by  direct  application  of  the 
elastic  solution  as  follows.  The  solution  to  the  elastic 
problem  with  the  coordinate  system  attached  to  the  crack 
tip, as  shown  in  Figure  12,  has  been  derived  by  Sih  (^2), 
where  the  stress  functions  introduced  in  sub-section  4.1  are 


V2j)  =  (?4c) 

J 

with  the  functions  g,(s  ,s2)  being 

J  1 

(1+s* ) 

g  (s  ,s  )  =  - 

1  1  2  tt[4siS2-(1+s2)2] 


tt(1+s2  )[4s  1s2-(l+s2  )2] 


(66) 


(67) 


(68) 


The  corresponding  displacement  functions  are,  through 
Equation  (63), 


1 3j ( s  j  »s2 )Q£n 


/i/z~Wc\ 
\  i  /zT  -  /c  / 


(69) 
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The  preceding  equations  represent  a  fundamental  solution  in 
that  they  may  be  used  to  construct  the  solution  for  an  arbi¬ 
trary  distributed  load  on  the  crack  surface.  For  the  problem 
of  interest  here,  it  is  necessary  to  find  the  solution  when 
uniform  tractions  of  intensity  p  are  applied  to  the  segment 
of  the  crack  lying  between  Xj  =  -d  and  Xj  =  0.  That  is,  let 


Q  -  pdC  and  c  =  -£ 


(70) 


and  integrate  Equation  (69)  over  the  interval  -d£^j<0,  the  result 
being 


W  = 


•p/ST  g-(s  , s  ) 
J  1  2 


( i  /FT  +  /cT 

2 /FT  -  —  (d+z  .)s,n  ^ - — 

J  /a  J  \  i  /FT  -  /cT 

*  J 


(71) 


The  stress  functions  $.(z.),  upon  differentiation  of  Eqs. 

J  J 

(71  )  are  found  to  be 


W  B  -pgj(s 


f 

) '  2/d/ z .  -  i £n 

J 


^lVzT  *  /d\| 

V  i  /z  •  -  /cT  / J 
J 


(72) 


The  solution  to  the  el  as t i c- pi ast i c  problem,  where  the 
loads  Q  are  at  a  distance  c  from  the  physical  crack  tip  and 
where  the  plastic  zone  extends  a  distance  d  beyond  the  crack 
tip  may  now  be  obtained  by  the  principle  of  superposition 
(see  Figure  14).  That  is,  replace  c  by  c  +  d  in  Equation  (66), 
let  p  =  -Y,  where  Y  is  the  yield  stress,  in  Equation  (72),  and 
the  stress  function  for  the  elastic-plastic  problem  is  equal 
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to  the  sum  of  Equations  (66)  and  (72).  The  displacement  function 
is  obtained  in  the  same  manner.  In  order  to  satisfy  the 
Dugdale  hypothesis  of  finite  stress  at  x:  =  0,  consider  the 


resultant  expression  for  4>.(z.)  which  is 

J  J 


l/c+S 


MZj)  -  -9.(3,  .s2){  - 2Y  S- 

3  3  3  /z-(z.+c+d)  Vzj 

V.  J  J 


i  Y£n 


i  /z  .  +  /? 

_ J 

i  /FT  -  /J 

J 


Thus,  it  follows  that,  for  non-singular  stresses  at  the  tip 
of  the  plastic  zone  (the  crack  tip  of  the  two  superposed 
elastic  solutions),  the  coefficient  of  1//FT  must  be  zero. 

J 

Then,  the  following  expression  relating  the  ratio  Q/Y  to 
the  geometric  quantities  must  hold: 


Q_  =  2- 

Yd  L 


The  expressions  for  <f>,(z-)  are  then 

J  J 


<f>j(Zj)  =  Yg^Sj.sJ  2/zTcT  -  i  ( d  +  z  j  )  £  n 


i /FT  +  /S 
_ J 

i  /zT  -  /cT 

J 


i /zT  +  /c+d 
_ J _ 

i  /FT  -  /c+d 
J 


This  displacement,  u„,  of  the  upper  crack  surface  may 

be  obtained  by  substitution  of  Equation  (75)  into  Equation  (62)  with 

z.  equal  to  x  .  The  result  is 
J  * 
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f 

* 


s  ( s 2  - 1  )  Y  • 
1  '  2  ' 


7TM  [4s  J  S  2  -  V 1  +S  |  )  ] 


<  2/- ( x  j d )  -  (d+x  x  )£n 


73  -  /^7 
/r  +  /-x , 


_  /c+d  -  /-x  \ 

+  2/d(c+d)£n  [  - 

,/c+Z  +  /^TT  / 


x<0 


(76) 


The  length  d  of  the  plastic  zone  is  independent  of 
velocity  and  is  related  to  Q/Y  and  the  distance  c  of  the 
load  Q  from  the  physical  crack  tip  by  Equation  (74).  The  re¬ 
lationship  is  shown  graphically  in  Figure  15. 

Since  the  plastic  zone  length  is  independent  of  velocity 
a  contour  C  which  surrounds  the  plastic  zone  and  is  fixed 
in  size  with  respect  to  time  may  be  used  in  calculating  the 
energy  dissipation  for  this  model.  Then,  the  path  indepen¬ 
dent  integral  around  the  contour  of  Equation  (8)  may  be  evalua¬ 
ted  by  shrinking  C  onto  the  boundary  of  the  plastic  zone. 

As  a  result,  the  expression  for  Gj  is 

dU  , 


'I 


3U. 

Ti  JT  ds  * 
1 


2  2  3X .  1 


(77) 


In  the  case  of  perfect  plasticity  considered  here, 
a22  is  equcl  to  Y,  the  constant  yield  stress.  Therefore 


au2 

r —  dx  =  Y [u  ]. 
3x  i  L  2jL 


(78) 


or 


G I  -  Y5T 


(79) 
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°a/o  onva  ssaais  cmi  01  avoi  asnddv 


Figure  15.  Finiteness  Condition  for  Semi  -  Inf i ni te  Crack. 
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LATIVE  PLASTIC  ZONE  LENGTH  d/c 


0* 

0 

*1 


where  6y  is  the  crack  opening  displacement  which  is  defined 
in  terms  of  the  displacements  u*  and  u^  of  the  upper  and 
lower  crack  surfaces  at  = -d  as 


6  =  [u*  -  u‘]  ■  (80) 

x  =  -  d 

i 

Then  making  use  of  Equation  (76),  the  expression  for 
for  the  semi -inf ini te  crack  model  is 


2  ( s  ^  - 1  )  Y  s  i 
Try[4sis2-(l+s^)2] 


r 

d  +  /d ( c+d )  £n 


/c+U  -  /cf  \ 
/c+d  +  /d  / 

(81) 


For  convi eni ence ,  using  the  above,  the  energy  release  rate 
may  be  written  as  the  product  of  the  normalized  elastic 
value  and  a  function  of  Q/YC.  That 


el 


f  (q) 


(82) 


with  q  =  Q/YC  and 

f(q)  =  ~~t  ^1  -q2  +  p£n  ( q  +  /I  +q2  )J  •  (83) 

Figure  16  is  a  plot  of  f(q)  versus  q.  Note  that  f(q)  is 

equal  to  the  ratio  of  Gj  to  the  corresponding  elastic  value 

Gt  ,  at  a  given  crack  speed. 

*el 
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NORMALIZED  LOAD  q  =  Q/Yc 


Figure  16.  Ratio  of  Plastic  -to  Elastic  Energy  Release 
Rate  as  a  Function  of  Normalized  Load  for 
Semi -Infinite  Crack. 
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4.4  CONSTANT  LENGTH  CRACK 


Consider  next  the  problem  of  a  crack  of  fixed  length 
2c  which  propagates  with  a  fixed  velocity  v.  This  crack 
model,  with  a  Dugdale  plastic  zone  is  shown  in  Figure  17. 
With  a=c, the  elastic  case  would  be  depicted. 


As  for  the  semi -i nf i ni te  crack  model,  the  general  re¬ 
sults  of  Appendix  I  apply  for  the  elastic  results.  Thus, 

Equation  (64)  also  applies  for  the  finite  crack  problem 
★ 

with  Gj  the  static  energy  release  rate  now  being 
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The  analysis  for  the  plastic  case  is  similar  to  that 

used  for  the  s emi - i nf i n i te  crack.  Using  the  results  of 

(12) 

Sih  solutions  for  the  sectionally  holomorphic  functions 

$.(z.),  for  a  crack  of  length  2a,  in  an  infinite  elastic 

J  J 

medium,  when  a12  =  0  for  all  x  on  the  crack  line,  are  of 
the  form 


W  = 


e'zj 


dC  + 


where  the  functions  f x ,  f  ,  T  and  ?2  are 

(1+S*) (at 2  +  U22 ) 

f.(0  =  — £ — - — 

(1+s2)2  -  4s  s 


4s  s 

f  U)  = - i-*- 

2  (1+s2) 


fxU) 


1 


°li  *  °22 


2(s=-s2) 

1  l 


2  2 


1  1 


+  a 


2  2 


1  +s 2 

2 


2K-s;> 


(88) 


(89) 


with  and  being  normal  stresses  on  the  upper  and 
lower  sides  of  the  crack  and  with  a”x  and  a*  being  the 
stresses  at  infinity. 

Consider  first  the  problem  where  concentrated  loads  Q 
are  acting  on  the  upper  and  lower  crack  surfaces  at  xl  =  +b, 
-b  (see  Figure  18).  For  this  case,  rx  =  T2  =  0  and 

2  ( 1  +  s 2 )  P 

f.(c)  = - [6(C+b)  +  6(5-b)].  (90) 

1  (1+s2 )2-4s  s 


Then  $.(z.)  are 

J  J 


tj(zj)  -  -2gj!s1,s2)Q 


/a2-b2  z  . 


(91) 


where 


S2)  are  given  by  Equations  (67)  and  (68). 


The  above  expression  may  be  used  to  generate  the  solu¬ 
tion  to  problems  where  uniform  tractions  are  imposed  upon 
the  crack  edge.  For  the  case  shown  in  Figure  19  where 
the  portions  of  the  crack  surface  extending  from  |xj  =  c 
to  | x  | =  a  are  under  pressure  p,  the  functions  $,(z.)  may 

A  J  J 

be  obtained  by  letting  Q  =  pd£  and  b  =  £  and  integrating 
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Figure  18. 


Constant  Length  Crack  Propagating  with 
Velocity  v  and  Loaded  by  Concentrated  Forces 


Figure  19. 


Constant  Length  Crack  Propagating  with 
Velocity  v  and  Loaded  Uniformly  over  Portion 
of  Crack  Surface. 


Equation  (91)  over  the  range  c^^a .  The  result  of  this  opera¬ 
tion  is 


W  ■  -Vs>-s*)p 


r 

2z ,cos" 1 (c/a) 

J 


/ z  . /az-cz  -  c /a2-z? 
-  i  log)  - - ■  ■  -  J 

\z  ./a2-c2  +  c/az-z? 
J  J 


+  7T 

J 


(92) 


In  order  to  obtain  the  solution  for  a  crack  length  2c 
with  a  plastic  zone  of  length  (a-c),  the  above  solution  with 
p  =  -v  is  added  to  the  elastic  solution  of  a  crack  of  length 
2a  under  uniform  tensile  stress  P  at  infinity  (see  Figure 
17).  The  latter  solution,  which  may  be  obtained  from 
Equations  (85)  through  (89)^^  is 


Vzj> 


=  9j(sI,s2)Pir 


H 


/zT=a*- 

.  J 


-  1 


+  aj 


(93) 


wi  th 


A 


(94) 


A 
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-i-  +  _ . 

1  +s |  2(Si-Sj) 


(95) 


Inspection  of  $.(z.)  for  the  two  cases  reveals  that, 

J  J 

for  the  superposition  of  the  two  solutions  to  satisfy  the 
Dugdale  hypothesis  of  finite  stress  in  the  neighborhood  of 
|xj=  a.  the  following  condition  must  be  met: 
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(96) 


The  superposed  solution  that  satisfies  this  condition  is 
then 


Vzj}  =  -gj(si-s*)Y  1 


1  og 


z . /a2-c2  -  c/a2-z?  \ 
— J _ I 

z ./a2-c2  +  c/a2-z^  J 

J  J 


(97) 


In  order  to  evaluate  di s p 1 acements, the  functions  <j).(z.) 

J  J 

are  obtained  bv  integrating  the  preceding  expression  with 
respect  to  z  . .  Thus  , 
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W  ■ 


-9j(VS2)Y  l2j’°9 


Zj/a2-c2  -  c/a2-z*  \ 


z  . /a 5 -c 2  +  c/a2 -z2. 
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(98) 


Then,  by  Equation  (62)  the  displacement  on  the  upper  por¬ 
tion  of  the  crack  surface  is 


u 


+ 

2 


S  !  (1  ~S2  )  Y 

ffp[4s  j  s2 - (1+s2  )2] 


j 1 og 
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(99) 
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The  crack  opening  displacement,  6-p  at  |x 


c  i  s 


4s  ( 1 -s2 ) Yc 

6t  = - - -  log(a/c) 

ttm[4s1s2-(1+s2)2] 


or,  making  use  of  Equation  (96) 


6 


T 


4s1(l-s2)Yc 
7TM  [4s  1s2-(l+s2)2] 


log[sec(|y-)]  . 


(100) 


(101) 


Use  of  the  Dugdale  hypothesis  has  led  to  a  relationship 
between  the  length  a  of  the  plastic  zone  plus  the  crack,  the 
length  c  of  the  crack,  the  applied  stress  P,and  the  yield 
stress  Y.  The  relationship  is  again  independent  of  crack 
speed  and ,  since  Equation  (9G)  is  the  same  as  Equation  (56),  is 
identical  to  the  static  result.  Thus,  calculating  the  energy 
release  rate  the  same  reasoning  as  was  used  for  the  case  of 
the  semi - i nf i ni te  crack  applies  and  the  general  expressions 
of  Equations  (77)  -  (79)  and  Equation  (82)  are  valid  in  this  case 
also,  with  f(q)  being  given  as, 

f (q )  =  — —  1  og ( sec [5-  q])  (102) 

7r 2  q  2  d 

where 

q  =  P/Y  .  (103) 

Again,  it  should  be  noted  that  f(q)  in  Equation  (102)  is  equal  to 
the  ratio  of  the  plastic  energy  dissipation  to  the  elastic 
energy  release  rate  for  the  Yoffg  model.  This  ratio. 
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independent  of  crack  speed,  is  plotted  in  Figure  21.  Thus, 

for  a  given  value  of  P/Y  the  value  of  GT  /G?  may  be  obtained 

*pl  1 

by  multiplying  the  plot  in  Figure  20  by  the  corresponding 
value  of  f(q)  obtained  from  Figure  21. 
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ENERGY  RELEASE 


PLASTIC  TO  ELASTIC  ENERGY  RELEASE  RATE 


SECTION  V 


PLASTIC  FLOW  AROUND  AN  EXPANDING  CRACK 

A  more  realistic  crack  model  for  the  projectile  pene¬ 
tration  problem  should  involve  a  central  crack  in  the  tar¬ 
get  material  moving  at  both  ends.  As  before,  the  material 
along  the  lines  of  expected  crack  extension  will  be  assumed 
to  flow  plastically  and  to  move  with  the  crack  but  at  a 
different  velocity.  The  energy  released  as  a  function  of 
the  crack  speed  will  be  examined. 

5.1  CRACK  EXPANDING  ELASTICALLY 

Broberg  considered  the  case  of  an  infinite  elastic 

plate,  under  uniform  tensile  stress  P,  with  a  crack  opening 
from  zero  length  and  extending  with  uniform  velocity  v. 

The  stresses  and  displacements  near  the  crack  tip  for  this 
case  are  given  in  Appendix  I.  It  should  be  noted  that  this 
model  would  apply  to  the  case  where  an  expanding  crack  has 
reached  a  length  several  times  greater  than  its  original 
length  prior  to  propagation  (after  the  initial  acceleration 
stage) . 

One  of  the  important  contributions  of  the  Broberg  solu¬ 
tion  was  that  it  allowed  a  more  realistic  estimate  of  the 
rate  at  which  energy  is  dissipated  at  the  crack  tip,  i . e .  , 

the  energy  release  rate  Gj.  Using  Equations  (14)  and  ( I -1 1 ) ,  a 

★ 

plot  of  the  normalized  energy  release  rate,  Gj/Gj,  versus 
the  ratio  of  crack  speed  to  dilatational  wave  velocity  is 
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given  for  steel  (v  =0.29)  in  Figure  26.  Two  factors  should 

be  kept  in  mind  when  using  this  plot.  First,  the  normaliza- 

* 

tion  with  respect  to  Gj  implies  that  the  value  of  Gj  at  a 
given  velocity  is  directly  proportional  to  the  instantaneous 
crack  length  c,  since,  with  the  value  kj  =  P/rr, 


Secondly,  the  value  determined  is  the  predicted  amount  of 
elastic  energy  irreversibly  absorbed  by  the  crack  tip  at  a 
given  velocity;  the  mechanism  by  which  this  energy  is  ab¬ 
sorbed  is  not  specified.  Thus,  in  order  to  derive  a  func¬ 
tional  relationship  between  the  material  parameters,  crack 
length,  load,  and  velocity,  this  mechanism  must  be  properly 
model ed . 

If,  for  example,  the  Griffith  criterion  of  direct  con¬ 
version  into  surface  energy  is  used,  the  velocity  could  be 
approximated  through  the  relation 

Gj  =  2y 

or 

Gj  1 6y y 

G*  ir(ic+l  )  P2c 

where  y  is  the  specific  surface  energy.  Thus,  the  predicted 
limiting  velocity  as  c-*»  would  be  the  Rayleigh  wave  velocity 
cR.  However,  the  above  criterion  can  only  be  applied  to  an 
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exceedingly  brittle  material  such  as  glass,  with  the  assump¬ 
tion  that  y  is  independent  of  crack  speed.  Experimental 
evidence  supports  the  fact  that  crack  velocities  in  ductile 
materials  are  dependent  upon  a  dissipative  process  and  may 
be  limited  to  speeds  which  are  only  a  fraction  of  c^  (29). 

5.2  PLASTIC  YIELDING  AROUND  THE  CRACK 

Incorporation  of  the  Dugdale  hypothesis  into  the  material 
of  dynamic  crack  propagation  allows  consideration  of  plastic 
energy  dissipation.  However,  as  will  be  seen,  this  does 
not  satisfactorily  specify  the  absorption  mechanism.  That 
is,  again  an  additional  condition  dealing  with  material  be¬ 
havior  is  necessary  if  one  desires  the  limiting  velocity 
for  a  gi^en  material  and  applied  stress.  The  information 
that  is  obtained  here  is  the  amount  of  energy  dissipated  by 
a  very  specialized  plastic  zone  for  a  given  crack  velocity. 

The  formulation  of  the  elastic  boundary  value  problem 

that  is  equivalent  to  the  elastic-plastic  problem  of  a 

uniformly  expanding  crack  with  a  Dugdale  plastic  zone  ahead 

of  the  crack  tip  is  made  in  the  same  fashion  as  for  the 

(14) 

steady-state  and  static  problems.  Atkinson  consi¬ 

dered  this  problem  for  a  crack  whose  tips  have  a  velocity 
v  with  a  plastic  zone,  the  end  of  which  is  traveling  at  ve¬ 
locity  8.  The  equivalent  elastic  problem,  then,  consists 
of  an  infinite  plate  under  uniform  tensile  stress  P  with  a 
crack  of  length  2a  =  2 6 1  upon  which  tensile  stress  Y  is 
applied  over  the  region  vt£|x,|£8t.  The  expressions  for 


normal  stress  °22.  ahead  of  the  plastic  zone  and  for  the 
second  derivative  with  respect  to  time  of  the  displacement, 
u2,  of  the  crack  surface  are 
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for  x2  =  0  and  |x2 |>60T^or  r>£2)  and 
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(105) 


where 


r  =  ctt  k  =  c2/cl 

b  =  v/ct  r  =  x2/t2  (106) 

60  -  0/c  i 

and  Y  is  the  yield  stress  of  the  material.  The  elastic 

wave  velocities,  c.  and  c  ,  are  described  in  Appendix  I. 

1  2 

The  functions  F(r)  and  tp ( r )  are 
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F(r)  = 


_ 4k4r32/l  -r _ 

( r- 2k 2 )  2-4k  V(  1  -r)  (  k2-i 
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for  B2<r<k2 
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(109) 


and  where  the  following  equation  may  be  used  to  obtain  B, 


+  m  + 1  j  ±1±  *(s)ds 

IT  IT  17  /=■ 


(110) 


An  algebraic  expression  for  B  in  terms  of  elliptic 
integrals  may  be  obtained  by  combination  of  Equations  (108),  (109), 
and  (HO).  The  result  obtained  by  Atkinson  has  been  cor¬ 
rected  here  and  is  as  follows 


BZi  =  2(P-Y)B0  -  ^  Z2 
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The  functions  K,  E,  and  n,  are  complete  elliptic  integrals 

of  the  first,  second,  and  third  ki nds, respecti vely,  and  are 

( 33  \ 

defined  in  Byrd  and  Friedman  v 

Equation  (105)  must  be  integrated  twice  in  order  to 
obtain  the  displacement  of  the  crack  line.  To  this  end. 
Equation  (105)  i :  rewritten  as 
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Application  of  Equation  (110)  simplifies  the  preceding  equa¬ 
tion  to  read 
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Recalling  that  the  expression  for  4>  ( s )  is  given  in 
Equation  (108),  the  integral  in  Equation  (115)  may  be  broken  down  into 
four  integrals  which  may  be  evaluated  by  well  known  methods 
of  complex  variable  theory  (30).  After  further  simplifi¬ 
cation  through  Equations  (107)  and  (109),  the  resulting  expression 

for  9 2 u  /9t2  is 
2 


9t2  pc2 \xi  j (Bj-r) 

0<r<B2 


Br- 


TAr- 


/Spr  (r-b2)/pr 


(116) 


For  any  function  g(r,x),  where  r  is  given  in  Equation  (106),  the 
chain  rule  of  elementary  calculus  allows  the  derivative  with 
respect  to  x  of  g  to  be  written 


9g  .  9g  9r  0  r ^2  9g 
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(117) 


71 


Thus,  Equation  (116)  may  be  written 
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Integration  of  the  above  expression  yields 
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Again  making  use  of  Equation  (117),  a  second  integration  with 
respect  to  r  leads  to  the  result 
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The  condition  that  stress  be  bounded  in  the  neighbor¬ 
hood  of  the  points  |x  |=B0t  is  satisfied  if  the  slope  of 
the  displacement  of  the  crack  surface  is  zero  at  |xi|=B  x. 
Differentiation  of  Equation  (120)  with  respect  to  x ^  reveals 
that  this  condition  requires  that 
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With  elimination  of  B  by  the  above  relation,  the  expression 
for  u2  may  be  written  in  terms  of  x  and  t  as 
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5.3  FINITENESS  CONDITION 

Although  Equation  (122)  does  not  explicity  depend  on  the 
applied  load  P,  there  is  implicit  dependence  through  the 
relationship  between  B  and  TA  of  Equation  (121).  This  equation 
together  with  Equations  (111)  through  (113)  imposes  a  functional 
dependence  between  the  parameters  P/Y,  Bq,  and  b.  Thus,  in 
evaluating  u,  only  two  of  these  may  be  independently  speci¬ 
fied.  Upon  combination  of  Equations  (111)  through  (113)  with  Equation 
(121),  the  finiteness  condition  may  be  represented  by  the 
equation 
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Using  the  precedi ng  equation,  P/Y  has  been  plotted  as 

a  function  of  relative  plastic  zone  size  (6  - b ) / 6  for  fixed 

o  0 

values  of  8Q  in  Figure  22  and  for  fixed  values  of  b  in 
Figure  23.  The  curve  for  80  =b=0  represents  both  the  static 

case  and  the  constant  velocity  fixed  length  crack.  An  im¬ 
portant  effect  to  be  noted  here  is  that,  for  a  given  value 
of  P/Y,  an  increase  in  crack  velocity  causes  a  decrease  in 
relative  plastic  zone  size.  Figure  24  shows  the  variation 
of  8 o ,  the  normalized  velocity  of  the  plastic  zone  tip, 
with  the  stress  ratio  P/Y  when  b  is  held  fixed.  It  can  be 
seen  here  that  8^  has  a  lower  limit  equal  to  the  normalized 
crack  speed  b  and  reaches  a  maximum  which  is  less  than  or 
equal  to  the  normalized  Rayleigh  wave  volocity  when  P/Y  =  l  . 

The  plastic  zone  tip  reaches  the  Rayleigh  wave  speed  only 
when  the  crack  tip  reaches  the  same  speed. 

5.4  ENERGY  RELEASE  RATE 

For  this  model,  the  rate  at  which  energy  is  dissipated 
by  the  moving  crack  tip  is  evaluated  by  the  result  of  Equation 
(20)  in  Section  II. 

Making  use  of  the  relation 
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Figure  22.  Finiteness  Condition  for  Uniformly  Expanding 

Crack.  Plot  of  Ratio  of  Applied  to  Yield  Stress 
versus  Relative  Plastic  Zone  Length  for  Various 
Velocities  of  the  Tip  of  the  Plastic  Strip. 
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CRACK  OPENING  DISPLACEMENT  tt^8t/|x|  Y 


Figure  25.  Plot  of  Normalized  Crack  Opening  Displacement 

Versus  Relative  Plastic  Zone  Length  for  Various 
Crack  Tip  Velocities.  Valid  for  Constant 
Length  Crack  and  Uniformly  Expanding  Model. 
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u  2 (x  j  >x )dx i  =  J  7— 1  dxi  +  u2(B0t,t)Bo  -  u2 (bx  ,x)b 

bx  / i o / 


(124) 


and  noting  that  u2(Bqx,x)  =  0,  Equatior  (20)  may  be  written  as 


G.b  =  2Yu  (bx,x)b  +  2Y  u  (x  ,x)dx 

1  2  3X  2  11 


(125) 


Upon  substitution  of  Equation  (122),  the  integral  in  the 
above  equation  may  be  expressed  in  the  form 
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(129) 


which  reduces  to  the  result 


I  =  t  2  b 2 
i 


[21og(^-)  - 


2/^TT^2 

_ _  0 
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' 't’J 


(130) 
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(x2-x2  b2 ) 1 og 


x/B^-b2  +  b/T232  -x2 
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or 


I  =  T2b/g2-b2  cos  (q— ) 

2  M0  P 

0 


(132) 


Therefore,  the  derivative  with  respect  to  x  of  Eq .  (126)  is 
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(133) 


From  Equation  (122),  the  expression  for  u2(bx,x)  is  obtained 
to  be 
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U2{bT,T) 


(bT)TA 

pc2{62-b2)V2b 
1  0 


log(|i) 


(134) 


and,  from  Equations  (125),  (133),  and  (134),  it  follows  that 
G j  is 


Gt  = 


2YTA ( bT ) 


r/iTrE‘  ,b 


pc2 ( B2 -b2 ) 3  2  b 
1  0 


cos"*  (H_)  -  log(^) 
o 


(135) 


The  energy  release  rate  normalized  with  respect  to  the 
corresponding  static  elastic  value  is  plotted  in  Figure  26 
and  may  be  compared  to  the  elastic  result  for  a  crack  of 
the  same  velocity  and  length. 


A  quantity  mentioned  previously  as  being  of  descriptive 
importance  is  the  crack  opening  displacement,  6y,  at  the 
point  |xj=vt  where  the  plastic  zone  begins.  This  quantity 
is  twice  the  displacement  uy(bt,b),  given  by  Equation  (134),  or 


6 


T 


2cTA _ 

pc2(82-b2)3'2b 
i  o 


1  o  g  ( 3  o  /  b )  . 


(136) 


Since  3  / b  =  a/c,  it  is  easily  seen  that  Equation  (136)  is  the 
same  as  Equation  (100),  the  crack  opening  diplacement  for  the 
fixed  length  case.  It  should  be  noted, however ,  that  for 
a  fixed  applied  stress  P,  the  oehavior  of  Sy  with  respect 
to  crack  velocity  is  not  the  same.  This  is  due  to  the  fact 
that  in  the  steady  state  case,  for  which  Equation  (100)  applies, 
the  ratio  a/c  is  dependent  only  on  P/Y  and  is  independent 
of  crack  velocity,  while  a  change  of  crack  velocity  in 
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Equation  (136)  when  P/Y  is  fixed  requires  a  change  in  B  /b  through 
the  finiteness  condition. 

The  relation  G  j  =  <5yY  has  been  used  in  Section  IV  for 
the  static  case.  Multiplication  of  Eq.  (136)  by  the  yield 
stress  Y  and  comparison  with  Eq.  (135)  reveals  that  this 
relationship  does  not  apply  here.  The  basic  difference  is 
that  expanding  crack  model  accounts  for  energy  dissipation 
not  only  through  the  translation  of  the  plastic  region  with 
the  crack  tip  but  also  through  enlargement  of  this  plastic 
region.  Therefore,  one  would  expect  that  Gj,  as  given  by 
Equation  (135),  would  be  higher  than  6jY.  This  is  shown  graphi¬ 
cally  ii  Figure  27  where  the  relative  difference  between 
Gj  and  6^-Y  is  plotted  as  a  function  of  normalized  crack 
velocity  b  and  in  Figure  28  where  the  same  difference  is 
plotted  as  a  function  of  P/Y. 
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CRACK  VELOCITY  v/c, 


Figure  27.  Relative  Difference  in  Energy  Release  Rate  for 
Expanding  Crack  Model  Due  to  Difference  in 
Method  of  Calculation  Plotted  as  a  Function  of 
Crack  Velocity. 
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Figure  28.  Relative  Difference  in  Energy  Release  Rate 
for  Expanding  Crack  Model  Due  to  Method  of 
Calculation  Plotted  as  a  Function  of  P/Y. 
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SECTION  VI 


! 

| 


I 


RESULTS 


General  results  have  been  presented  in  Section  II  thar 
allow  the  calculation  of  the  rate  at  which  energy  is  released 
by  a  uniformly  moving  crack.  Of  particular  interest  is  the 
fact  that  a  path  independent  integral  may  be  used  for  this 
calculation  not  only  for  the  steady  state  cases  of  crack 
propagation  but  for  the  non-steady  state  case  of  the  uni¬ 
formly  expanding  crack,  as  modeled  by  Broberg. 


The  rate  at  which  energy  is  dissipated  due  to  plasti¬ 
city  was  calculated,  for  the  dynamic  case,  in  a  manner 
analagous  to  the  analysis  of  Goodier  and  Field  for 

the  static  case.  Although  a  specific  result  was  presented 
only  for  perfect  plasticity,  work  hardening  can  be  incor- 
pora*-'1  through  Equation  (17). 


',rder  to  gain  further  insight  into  the  transient 
response  of  the  cracked  plate, it  is  useful  to  review  the 
wave  patterns.  When  a  transient  load  is  applied  to  the 
crack  surface,  the  result,  aside  from  wave  propagation  from 
the  loaded  surface, is  that  the  crack  tips  form  the  center 
of  two  outgoing  cylindrical  waves.  For  the  period  of  time, 
t  ,  before  these  waves  begin  to  interact  the  short  time  be¬ 
havior  of  the  solution  is  described  by  the  equations 


a 


i  j 


k  ) 

-  fij{9) 

SZr  1J 


0<t£t$ 
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where  r  is  the  distance  from  the  crack  tip,  and 


I 


{ 


i 

i 


The  function  f..(6)  is  the  same  as  in  the  static  case.  Since 
i  J 

the  disturbance  at  one  crack  tip  is  not  yet  affected  by  the 
presence  of  the  other  crack  tip,  this  is  the  range  for  which 
the  solution  to  the  problem  of  sudden  loading  on  a  semi¬ 
infinite  crack  is  valid.  For  normal  loading  the  stress- 
intensity  factor  may  be  obtained  from  the  solution  of  Baker 
^ '  '  which  is, 

k|s)  =  o  / a  . 785/c  t/a) 

I  1  2 

for  the  case  of  steel  . 

This  result  is  plotted  in  Figure  9  and  shows  good 
agreement  with  the  numerical  solution  obtained  here. 

The  stress  wave  pattern  loses  its  geometrical  simpli¬ 
city  after  the  two  cylindrical  waves  begin  to  interact  al¬ 
though,  very  close  to  the  crack  tip,  the  angular  distribu¬ 
tion  and  square  root  singularity  remain  the  same.  It  is 
during  this  period  of  maximum  disorder  that  the  stress- 
intensity  factor  takes  on  a  maximum  value.  The  pattern 
described  here  is  shown  in  Figure  7  for  various  times 
after  the  initial  loading.  The  outer  wave  front  surrounding 
the  crack  is  not  to  scale  and  is  obtained  by  allowing  a 
sufficient  time  to  elapse  such  that  c2t>>2a.  At  this  time 
enough  energy  has  been  transmitted  away  from  the  crack  that 
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the  dynamic-stress  intensity  factor  is  very  close  to  the 
static  value. 

These  results  presented  here  are  in  agreement  with 
those  of  Soltesz  and  Sommer^^who  experimentally  deter¬ 
mined  the  time-dependent  stress-intensity  factor  due  to 
normal  loading  of  a  crack  in  a  finite  plate  of  Araldit  B. 
Measurements  of  the  crack  opening  displacement,  which  is 

directly  related  to  the  stress- intensity  factor,  showed 

_  5 

that  ki  reached  a  maximum  at  approximately  10  seconds, 
the  same  order  of  magnitude  as  obtained  here.  The  oscil¬ 
lations  that  were  observed  in  crack-opening  displacement 
were  explained  as  due  to  a  combination  of  wave  reflection 
at  the  edges  of  the  finite  plate  and  oscillation  in  the  k 
value  due  to  the  initially  applied  load,  i.e.,  the  result 
demonstrated  in  this  work. 

The  results  of  Section  IV  for  the  steady-state  cases 
may  present  some  difficulty  when  an  attempt  is  made  to  in¬ 
terpret  them  physically.  This  may  be  seen  in  the  elastic 
case  by  examination  of  Figure  20  which  applies  for  both 
the  semi - i nf i ni te  and  constant  length  crack.  If  the  assump 
tion  is  made  that  the  crack  will  run  at  a  speed  that  causes 

Gj  to  be  equal  to  a  critical  value  Gjc>  then, for  a  given 

★ 

velocity, an  increase  in  Gj  ,  the  corresponding  static 
value  of  Gj,  would  have  the  effect  of  lowering  the  crack 
velocity  in  order  to  maintain  the  same  critical  value  of  Gj 
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In  the  case  of  the  s emi - i nf i ni te  crack, it  can  be  argued 

that  the  result  of  an  increase  in  the  applied  load  P  would 

be  to  increase  the  distance  c  from  the  crack  tip,  thus  pos- 

* 

sibly  maintaining  the  same  value  of  Gj.  The  results  for 
both  the  elastic  and  plastic  case  should  have  some  use  in 
considering  the  physical  problem  of  a  wedge  load  applied  to 
the  crack  surface. 

The  character  of  the  elastic  stress-field  ahead  of  a 
moving  crack  tip  is  the  same  for  all  constant  velocity  models 
of  the  same  mode  as  stated  in  Appendix  I.  Thus,  since  the 
Yoffe  model  was  the  first  dynamic  crack  model, the  results, 
in  this  respect,  were  a  valuable  contribution.  However,  it 
is  important  to  note  that  the  model  of  a  crack  which  opens 
at  one  end  and  closes  at  the  other  is  physically  unrealizable, 
and  one  should  expect  paradoxes  in  the  result  of  any  analysis. 
The  fact  that  the  results  indicate  an  increase  in  the  elastic 
energy  release  rate  with  velocity  for  a  given  load  and 
crack  length  is  such  a  paradox  since,  in  effect,  it  says 
that  for  a  given  load  the  larger  the  crack  length  the  lower 
the  corresponding  crack  speed.  Further,  note  in  the  elastic 
case  that  the  energy  release  rate  for  the  entire  system  is 
zero  (no  energy  is  transferred  over  a  boundary  surrounding 
and  moving  with  the  entire  crack)  while  the  energy  release 
rate  for  one  crack  tip  approaches  infinity  as  the  crack 
speed  approaches  the  Rayleigh  wave  velocity.  In  the  corres¬ 
ponding  Dugdale  model,  the  length  of  the  plastic  zone  is 
independent  of  the  crack  velocity.  This,  in  itself  does 


not  agree  with  the  analysis  of  the  more  realistic  Broberg 
model.  In  addition.  Equation  (99)  indicates  that  the  displace¬ 
ments  vary  from  zero  to  infinity  along  the  crack  surface 
when  the  crack  speed  approaches  c^.  This  violates  the  as¬ 
sumptions  of  infinitesimal  strain  theory  and  indicates  the 
solution  is  invalid  for  large  values  of  crack  velocity. 

It  is  the  conclusion  here  that  the  Yoffe  model  should 

not  be  used  to  evaluate  the  quantitative  effect  of  crack 

velocity  on  energy  dissipation  or  plasticity  and  strain 

/  o  9 \ 

rate  effects.  Kanninen  et  al  v  '  used  this  model  in  deal¬ 
ing  with  strain  hardening  and  strain  rate  effects.  Their 
results  are  valid  for  low  velocities  since  the  value  used 
for  crack  surface  displacement  is  dependent  on  velocity 
through  the  function  plotted  in  Figure  20.  These  authors, 
in  fact, only  used  the  Yoffe  model  as  supporting  evidence 
for  the  approximation  of  the  dynamic  case  by  the  static 
values.  Inspection  of  the  corresponding  curves  for  the 
Broberg  model  i nd i cates  that ,  while  there  is  a  range  over 
which  the  quasi-static  case  may  be  used,  it  is  not  as  large 
as  anticipated  by  the  results  of  Figure  20. 

The  Broberg  model  corrects  many  of  the  defects  found 

in  the  Yoffe  model.  It  more  closely  approximates  the 

physical  situation  in  u.at  it  takes  into  account  the  fact 

the  crack  is  increasing  in  size.  The  elastic  and  plastic 

results  for  the  strain  energy  release  rate  (Figure  26) 

★ 

follow  a  trend  opposite  to  the  Yoffe  model,  with  Gj/Gj 
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approaching  zero  as  v+c^.  Therefore,  although  the  results 
for  both  models  have  indicated  the  Rayleigh  wave  speed  as 
an  upper  limit  on  crack  velocity,  the  reasons  are  exactly 
opposite.  Physically,  tne  arguments  for  a  limiting  velocity 
based  on  either  creation  of  surface  energy  (see  sub-section  5-1) 
or  plastic  dissipation  make  more  sense  when  used  in  conjunc¬ 
tion  with  the  Broberg  model. 


When  the  Dugdale  criterion  is  incorporated  into  the 

Broberg  model,  the  finiteness  condition  dictates  that,  for  a 

fixed  value  of  P/Y,  larger  normalized  crack  velocities  b 

decrease  the  relative  plastic  zone  size  (B  -b)/B  (see  Fig- 

0  0 

ure  23).  In  addition,  at  high  velocities  the  plastic 
value  of  energy  dissipation  approaches  the  elastic  value 
(Figure  26).  These  results  are  supported  by  experimental 
evidence  which  indicates  that  high  velocity  fracture  is 
essentially  a  brittle  process. 
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APPENDIX  I 


ELASTIC  PROBLEMS  OF  UNIFORM  CRACK  EXTENSION 


In  terms  of  polar  coordinates  r  and  0  with  the  origin 
at  the  crack  tip  as  shown  in  Figure  1,  Sih  ^  has  shown 
that,  for  a  class  of  dynamic  crack  propagation  problems  in 
which  the  crack  tip  moves  at  a  constant  velocity  v,  asymp¬ 
totic  expansions  of  the  dynamic  stresses  and  displacements 
as  r+0  can  be  found  without  specifying  crack  geometry  and 
loading  conditions. 

For  the  opening  mode  of  crack  extension,  where  the  loads 
are  placed  symmetrically  with  respect  to  the  crack  line,  as 
in  the  Yoff£  and  Broberg  models,  the  singular  portion  of 
the  dynamic  stress  field  near  the  crack  tip  is 


a 

1 1 


a 


2  2 


kg 

=  -  F  (s  ,s)[(l+s|)(2s2+l-s2)f(s  ) -4s  s  f(s  )] 

/27  1  1  2  121122 

(I'D 

—  F  (s  , s  ) [4s  s  f(s  )  -  (l+s2)2f(s2)]  (1-2) 

✓Zr  1  1  2  1  2  2  2  1 


12  m 


2s  ( 1 +s 2 ) F  (s  , s  )[g(s  )-g(s  )]  . 
1  2  112  1  2 


(1-3) 


The  corresponding  displacement  expressions  are 


uu  =  k  /2r  F  (s  ,s  ) [ ( 1 +s 2 ) ( cos 2 6  +  s2sin20)f(s  ) 
11112  2  1  1 

+  2s  s  (cos20  +  s2sin20)f(s  )] 

12  2  2 
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(1-4) 


PU2  =  ki/2r  F J ( s J ,s2)[si(1+s^)(cos20  +  s 2 s i n2 9 )g ( s j ) 

-  2si(cos20  +  s2sin28)g(s2 )]  (1-5) 

where  kj  is  the  opening  mode  s tress - i ntensi ty  factor  for  the 

corresponding  static  case.  The  functions  s.  (j  =  l  ,2)  are 

J 

Sj  =  [1  -  (v/Cj)2]^  j  =  1 ,2  (1-6) 

in  which  ct  and  c2  are, respecti vely, the  dilatational  and 
shear  wave  velocities  in  an  infinite  elastic  continuum.  In 
terms  of  the  shear  modulus  p,  Poisson's  ratio  v,  and  the 
mass  density  p,  Cj  and  c2  are 

[2y(l - v ) / p ( 1 -2v)]^  ,  plane  strain 

ci  =  ' 

[2p/p( 1 -v) ,  plane  stress 
c2  =  (p/p)^  • 

The  quantities  f(s.)  and  g(s.)  describe  the  angular  distri- 

J  J 

bution  of  the  stresses  and  displacements  and  are  obtainable 
from  : 

f2(s  ■)  +  g 2  ( s  . )  =  sec0(l+sUan20)'^ 

J  J  J 

( j  =  1 » 2 )  (1-9) 

f2(s.)  -  g 2  ( s  . )  ~  sec0(l  +sHan26)"  1 

J  J  J 


(1-7) 

(1-8) 
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The  function  F1(s1s2)  depends  on  the  particular  problem 
under  consideration.  For  the  crack  of  fixed  length  propa¬ 
gating  at  a  uniform  velocity,  FjCs  ,s  )  1S 

F  x  ( s  ! ,  s  2 )  =  l^s^  -  (1+s2)2]'1  .  ( I  - 1  0  ) 


For  Broberg's  solution. 


F  (s  ,s  )  =  s  /[(1+s2,-  -  4s2s2]K(s  )  -  4s2(1-s2)K(s  ) 
112  1  1  2  12  1  1  2  2 


[ 4 s 2  +  ( 1 +s 2 ) 2 ] E ( s  )+8s2 


?E(s2)}-. 


(1-11) 


The  expression  for  k^  is  seemingly  the  same  for  both 
models,  that  is 


k  =  p/c 
i 

where  c  is  one-half  the  crack  length.  However,  for  the  Yoffd 
model  c  is  constant  whereas  for  the  Borberg  model  c  is  a 
function  of  t.  Thus,  in  calculating  energy  release  rate, 
one  should  be  aware  of  the  restrictions  in  sub-section  2.3. 

For  a  semi -i nf i ni te  crack  with  the  load  and  geometric 
configuration  of  Figure  12,  F  (s,,sj  is  also  given  by 
Equation  ( I -1 0)  and  the  stress  intensity  factor  k  is  equal  to 
[(Q/tt)/(2/c)]  . 
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APPENDIX  II 
NUMERICAL  PROCEDURE 


The  dual  integral  equations  formulated  in  Section  III 
were  reduced  to  a  Fredholm  integral  equation  that  must  be 
solved  numerically.  This  integral  equation  is  in  the  form 


A*U,<)-f  A*(e,K)K{£,K)dn  =  h(S) 


where  h(£)  and  K(£,n)  are  known  functions. 


(II-D 


The  method  of  solution  that  was  chosen  was  to  approximate 
the  integral  in  Equation  ( 1 1 -1 )  by  a  finite  sum  and  solve  the  fol¬ 
lowing  set  of  simultaneous  linear  equations  for  the  unknowns 
A  ( £  i ,  < ) . 


A  Uis<)-  J  a  j  K(Ci ,  n  j )  A  ( n  j  ,<)  =  hU.) 

j —  1 


(II-2) 


i  =  !,••*,  N  . 


The  constants,  a.,  are  chosen  in  accordance  with 
Simpson's  rule  and 


£  .(ki> 

tji  { N- 1  ) 


nj  s  -HttI  • 


( 1 1-3) 


The  function  A  is  thus  approximated  for  each  choice 
of  the  parameter, k  ,  which  is  related  to  the  Laplace  trans¬ 
form  variable  p  by  the  equation, 


( 1 1 -4  ) 
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By  this  means,  for  a  given  value  of  K,  the  function  A  (£,k) 

can  be  represented  numerically  in  the  Laplace  transformed 

★ 

domain.  The  curve  showing  A  (1,k)  as  a  function  of  k  is 
shown  in  Figure  8. 


In  order  to  complete  the  solution  for  the  near-field 
stress  components,  the  function  M(t)  defined  by  Equation  (53) 
must  be  evaluated.  Noting  that 


L 


t 

r 


m(t)dt 


(1 1-5) 


it  follows  from  Equations  (53)  and  (51)  that 

M(t)  =  L~ 1  [f*(p)A*(l ,<)].  (II-6) 


For  the  particular  example  treated  numeri cal ly, f (t )  is  the 

'ft 

Heaviside  step  function  and  therefore  f  (p)  is  1/p.  The 

solution  is  then  completed  by  evaluating  numerically  the 

★ 

inverse  Laplace  transform  of  the  function  A  (l,c2/pa)/p. 
The  Laplace  transform  pair  in  question  may  now  be  written 


★ 

A  (1 ,c2/pa) 
P~ 


00 

► 

M(t)exp(-pt)dt 

0 


M(t) 


★ 

A  (1  ,c2p/a) 

- - -  exp(pt)dp  . 


(II-/) 


( 1 1-8 ) 


Make  the  following  definition 


c  =  PA  =  1 

C  K 
2 


( 1 1-9 ) 
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Ajd.s) 


=  A  (1,1/S) 


M  (T)  =  M(t)  =  M(aT/c  ) 
1  2 


(11-10) 


(11-11 ) 


(11-12) 


Then,  in  terms  of  the  variables  s  and  T,  the  Laplace 
transform  pair  is 


A* ( 1  ,  s 


l«  /t\ 

n  l  ( T )  = 


M  (T)exp(-sT)dT 

i 

6 

-  j  --1 (y?—  exp ( sT )ds 


(11-13) 


(11-14) 


Having  available  the  solution  of  Equation  (II-D,  A  (1  ,s)/s  is 

known  for  discrete  values  of  s.  There  are  a  number  of 

methods  available  that  allow  calculation  of  a  function  from 

knowledge  of  a  finite  number  of  values  of  its  Laplace  trans- 

1 35 ) 

form.  The  one  used  here  is  due  to  Papoulis'  as  en¬ 


larged  upon  by  Miller  and  Guy 


In  this  case , 


A  (l,s)/s  is  evaluated  at  points  given  by 


s  =  (l+n)S  ,  n  =  0,1  ,2,  *  *  * . 


(11-15) 


where  <5  i  s  a  real  positive  number.  In  order  to  change  the 
infinite  range  of  integration  in  Equation  (11-13)  to  a  finite  one, 
the  following  definitions  are  made: 


i 


•’  wax  „ 


x  =  2exp(-<5T)  -  1 


i|>(x)  =  M,(T)  =  Mx(  -  \  An  [  ( x  +  1  )/2]  ) 


Equation  (11-13)  is  then  changed  to 


A*(l  ,  [1 +n]6 )  i  | 


(1+n) 


,n  +  l 


( 1  +x )  %  ( x )  dx 


The  function  <Mx)  is  expanded  in  a  series  of 
nomials  orthogonal  on  the  interval  (-1,1)  as 


'l'(x)  =  l  C.P.  (x) 
k=0  K  K 


where 


p„(x>  =  ^ -  [(i-x)nn+x)n] 


n  2 n n  !  dx11 


and , 


+  1 


P„(x)P(x)dx 
n  m 


n^m 


n=m 


v  n+2 


Note  that 


0+x)n  =  l  a,P.(x) 
k  =  0 


(11-16) 

(11-17) 


(11-18) 

Legendre  poly- 
fol 1 ows 

(11-19) 


(11-20) 


(11-21  ) 


(11-22) 


where 
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(11-23) 


,  _  on / OL-  +  1  \  n(n-l  )  •  • « [n- ( k- 1  )] 
k  ‘  1  }  ( n+1 ) ( n+2 ) • • • ( n  +  k+1  ) 


Then,  making  use  of  Equations  (11-21)  through  (11-23),  sub¬ 
stitution  of  Equation  (11-19)  into  Equation  (11-18)  yields  the  follow¬ 
ing  set  of  equations  that  may  be  solved  for  C^: 


A* ( 1 ,[l+n]5) 
( 1+n ) 


?  n ( n- 1 )  ♦  • » [n-  (  k- 1  )  ]  Q  (n_24) 
k  =  0  (n+1 ) (n+2) • • • (n+k+1  )  k 

n  =  l ,2  ,  •  •  • , 


Note  that,  for  n=0, 


A  *  ( 1  , 6  )  =  C  • 

i  o 


(H-25) 


If  N  coefficients  are  calculated, a  partial  sum  is  ob¬ 
tained  from  Equation  (11-19),  and  by  Equations  (11-16)  and  (11-17)  an 
approximation  of  M^(T)  is  given  as 

N-l 

Mj(T)  =  l  C.Pk[2exp(-6T)-l]  .  (11-26) 

k  =  0 

The  parameter  6  is  chosen  so  that  M(T)  is  best  described 
for  a  particular  range  of  T. 
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12.  SPONSORING  Mill  T  ARY  ACTIVITY 


Air  Force  Armament  Laboratory 
Air  Force  Systems  Command 


This  report  is  concerned  with  the  application  of  fracture  mechanics  theory  to 
the  problem  of  projectile  penetration  and  the  situation  where  an  excessive  amount 
of  energy  is  available  to  damage  the  target  material  in  the  form  of  fracture  followed 
by  crack  propagation.  The  conditions  under  which  this  phenomenom  takes  place  depend 
on  several  parameters,  such  as  the  speed  and  shape  of  the  projectile  and  the  geometry 
and  material  of  the  target.  These  conditions  have  been  simulated  in  the  laboratory. 
The  present  investigation  deals  only  with  the  development  of  an  analytical  fracture 
model  of  the  projectile  penetration  problem  using  the  theory  of  continuum  mechanics. 

The  amount  of  energy  required  to  produce  a  crack  of  certain  size  in  a  pre¬ 
stretched  plate  within  a  very  short  period  of  time  is  estimated  by  solving  a  system 
of  dual  integral  equations.  Computer  programs  are  being  developed.  In  addition, 
the  energy  supplied  to  maintain  a  running  crack  at  a  certain  speed  is  obtained  with 
results  illustrated  by  graphs.  The  effects  of  plasticity  and  rate  sensitivity  of 
the  material  are  also  discussed  in  this  report. 
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